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Abstract. Spherically symmetric (ID) and two-dimensional (2D) supernova simulations for progenitor stars between 11 M© 
and 25 M© are presented, making use of the PrometheusA^ertex neutrino-hydrodynamics code, which employs a full spectral 
treatment of neutrino transport and neutrino-matter interactions with a variable Eddington factor closure of the (9(v/c) moments 
equations of neutrino number, energy, and momentum. Multi-dimensional transport aspects are treated by the “ray-by-ray plus” 
approximation described in Paper I. We discuss in detail the variation of the supernova evolution with the progenitor models, 
including one calculation for a 15 M© progenitor whose iron core is assumed to rotate rigidly with an angular frequency of 
0.5 rads“^ before collapse. We also test the sensitivity of our 2D calculations to the angular grid resolution, the lateral wedge 
size of the computational domain, and to the perturbations which seed convective instabilities in the post-bounce core. In 
particular, we do not find any important differences depending on whether random perturbations are included already during 
core collapse or whether such perturbations are imposed on a ID collapse model shortly after core bounce. Convection below 
the neutrinosphere sets in 30-40 ms after bounce at a density well above 10^^ g cm“^ in all 2D models, and encompasses a layer 
of growing mass as time goes on. It leads to a more extended proto-neutron star structure with reduced mean energies of the 
radiated neutrinos, but accelerated lepton number and energy loss and significantly higher muon and tau neutrino luminosities 
at times later than about 100 ms after bounce. While convection inside the nascent neutron star turns out to be insensitive 
to our variations of the angular cell and grid size, the convective activity in the neutrino-heated postshock layer gains more 
strength in better resolved models. We find that low (I = 1,2) convective modes, which require the use of a full 180 degree 
grid and are excluded in simulations with smaller angular wedges, can qualitatively change the evolution of a model. In case 
of an 11.2 M© star, the lowest-mass progenitor we investigate, a probably rather weak explosion by the convectively supported 
neutrino-heating mechanism develops after about 150 ms post-bounce evolution in a 2D simulation with 180 degrees, whereas 
the same model with 90 degree wedge fails to explode like all other models. This sensitivity demonstrates the proximity of our 
2D calculations to the borderline between success and failure, and stresses the need to strive for simulations in 3D, ultimately 
without the constraints connected with the axis singularity of a polar coordinate grid. 
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1. Introduction 

The mechanism by which massive stars explode is still unclear. 
State-of-the-art models with a spectral treatment of the neu¬ 
trino transport by solving the Boltzmann equation or/and its 
moments equations agree in the finding that in spherical sym¬ 
metry (ID) neither the prompt bounce-shock mechanism nor 
the delayed neutrino-driven mechanism lead to explosions for 
progenitors more massive than about lOM© (e.g., Rampp & 
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Janka 2002; Liebendorfer et al. 2001, 2004; Thompson et al. 
2003; Sumiyoshi et al. 2005). Previous two-dimensional (2D) 
simulations (e.g., Herant et al. 1994, Burrows et al. 1995, 
Janka & Muller 1996, Fryer 1999, Fryer & Heger 2000) and 
three-dimensional (3D) models (Fryer & Warren 2002, 2004) 
show the importance of convective overturn in the neutrino¬ 
heating layer behind the stalled supernova shock, which can 
enhance the energy transfer from neutrinos to the stellar matter 
and thus cause “convectively supported neutrino-driven explo¬ 
sions”. These multi-D models, however, employed radical sim¬ 
plifications of the treatment of neutrinos, mostly by grey diffu¬ 
sion or in a parametric way as heating terms. Concerns about 
the reliability of such approximations of crucial physics in stud- 
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ies of the supernova explosion mechanism were expressed by 
Mezzacappa et al. (1998). 

Also the influence of convective activity inside the nascent 
neutron star, i.e. below the neutrino sphere, on the explosion 
mechanism has long been a matter of debate and requires fur¬ 
ther studies. The Livermore group (Wilson & Mayle 1988, 
1993) obtained explosions in their basically ID models by as¬ 
suming that so-called neutron-flnger mixing instabilities exist 
in the newly formed neutron star, which accelerate the energy 
transport from the neutron star interior to the neutrino sphere. 
Thus the neutrino luminosities are boosted and the neutrino 
heating behind the supernova shock is enhanced. The analy¬ 
sis by Bruenn & Dineva (1996) and more recently by Bruenn 
et al. (2004), however, has demonstrated that neutrino diffusion 
leads to lepton number equilibration between perturbed fluid 
elements and their surroundings that is faster than assumed by 
Wilson & Mayle (1988, 1993). Therefore neutron Angers are 
unlikely to occur in the supernova core. Bruenn et al. (2004) 
instead discovered a new mode of doubly-diffusive instability, 
which they termed “lepto-entropy Angers” and which is also as¬ 
sociated with neutrino-mediated thermal and lepton diffusion. 
The importance of this phenomenon during the early, critical 
phases of the explosion, however, was recently questioned by 
Dessart et al. (2005) because of its slow growth compared to 
Ledoux convection. The latter, in turn, was predicted to play 
a role in supernovae on grounds of ID models of the neutrino 
cooling phase of nascent neutron stars. A Ledoux-type of con¬ 
vection was indeed found to be present during the first second 
after neutron star formation in 2D hydrodynamic simulations 
by Keil (1997), Keil et al. (1996), Janka & Keil (1998) and 
Janka et al. (2001). The latter simulations, however, considered 
only the proto-neutron star without self-consistently follow¬ 
ing its feedback with the environment of the supernova core. 
Moreover, a grey, flux-limited equilibrium “ray-by-ray” diffu¬ 
sion code for the neutrino transport was used, with strong sim¬ 
plifications in the description of the opacities. 

Only recently multi-dimensional simulations of stellar core 
collapse and post-bounce evolution with a spectral treatment 
of the neutrino transport have become possible (Buras et al. 
2003, 2006; Livne et al. 2004; Walder et al. 2005; Swesty & 
Myra 2005a,b; Burrows et al. 2006). Although these current 
approaches are the first steps of removing the severe deficien¬ 
cies of the previous generation of multi-dimensional models, 
all of them still contain approximations of various, and differ¬ 
ent, aspects in the treatment of 2D transport. Swesty & Myra 
(2005a,b), for example, use a flux-limited diffusion descrip¬ 
tion, an approximation also made by Walder et al. (2005) and 
Dessart et al. (2005), who in addition solve the transport for 
all neutrino energy groups independently. In contrast, Buras 
et al. (2006) have developed a “ray-by-ray plus” approximation 
based on a variable Eddington factor solver for the coupled set 
of neutrino moments equations and Boltzmann equation, in¬ 
cluding a full coupling of the energy bins by neutrino reactions 
and by Doppler and gravitational redshift effects. 

The approximations employed by the different groups are 
diverse and might hamper a detailed quantitative comparison of 
the results in the near future, and might constrain such efforts 
to a purely qualitative level. Eventually it will be necessary to 


test and possibly replace the current approximations by a more 
rigorous solution of the transport problem in the five- or six¬ 
dimensional phase space and in a relativistic framework, once 
the corresponding codes have become available and the nec¬ 
essary substantial increase of computer power has happened 
(Cardall 2004, Cardall et al. 2005). 

Here we present results obtained with the multi¬ 
dimensional neutrino-hydrodynamics code MuDBaTH, 
which is the “ray-by-ray plus” implementation of the 
Prometheus/Vertex code described in detail in Buras et al. 
(2005; Paper I). In continuation of our previous work (Buras 
et al. 2003, 2006), where also a broader introduction into 
the status of the field and its open questions is provided, we 
present here ID simulations for nine different progenitor stars 
with masses between 11.2Mo and 25 M©, and compare them 
with 2D simulations for three of these stars. The core collapse 
and post-bounce evolution of these models was followed until 
nearly 300 ms after shock formation. 

Using a state-of-the-art treatment of spectral neutrino trans¬ 
port for hydrodynamical supernova simulations, the main goals 
of our work are: 

- We compare ID and 2D models in order to obtain quanti¬ 
tative information about the influence of convection below 
the neutrino sphere on the neutrino emission, the evolution, 
and the structure of the nascent neutron star. We analyse the 
influence of proto-neutron star convection on the conditions 
in the neutrino-heating layer behind the shock, and assess 
quantitatively the impact of convective activity in the post¬ 
shock layer on the possibility for reviving the stalled shock 
and for getting a delayed explosion. 

- We also study the differences of convection between a non¬ 
rotating and a rotating 15 M© model, whose iron core spins 
rigidly before collapse with a period of about 12 seconds, 
leading to an “extreme” period of the settled neutron star of 
the order of 1 ms. In addition, we investigate the effects of 
low-mode (dipolar, I = 1, or quadrupolar, I = 2) hydrody¬ 
namic instabilities during the post-bounce evolution of an 
11.2M© star, comparing simulations with a full 180° polar 
grid and simulations which are contrained to a ~ 90° equa¬ 
torial wedge (with periodic angular boundary conditions), 
thus preventing the development of such low modes in the 
pattern of the fluid flow. 

- Moreover, we perform tests for the influence of (i) the nu¬ 
merical resolution, in particular in the lateral direction of 
our 2D polar grid, (ii) of the chosen size of the angular 
wedge, and (iii) of the way in which we perturb our mod¬ 
els to initiate the growth of convective instabilities, i.e., 
whether we follow a perturbed 2D model through core col¬ 
lapse and core bounce, or whether we map a ID model to 
the 2D grid shortly after bounce, imposing random pertur¬ 
bations at that time. 

The paper is structured in the following way. Main re¬ 
sults of our ID supernova simulations for the chosen set of 
progenitor models — whose basic properties are compared in 
Appendix A — will be discussed in Sect. 2, supplemented with 
more details in Appendix B. The 2D models will be presented 
in Sect. 3, with an analysis of the effects of convection in the 
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forming neutron star in Sect. 3.1, a discussion of convection 
in the neutrino-heating layer in Sect. 3.2, a description of our 
full 180° model in Sect. 3.3, of the rotating model in Sect. 3.4, 
and of neutrino emission anisotropies in Sect. 3.5. Most of 
our 2D simulations without rotation were started from ID col¬ 
lapse models only shortly after bounce, at which time small 
random perturbations were imposed to seed convective insta¬ 
bilities. Since the adequacy of such an approach may be dis¬ 
puted, we also performed simulations where the collapse phase 
was followed in two dimensions. This allowed us to investi¬ 
gate the growth of inhomogeneities during infall and to assess 
the possible influence of that on the growth of convection after 
bounce (Appendix E). A summary and conclusions will follow 
in Sect. 4. Appendix C contains a linear analysis of the struc¬ 
tural changes of the proto-neutron star which can be expected 
as a consequence of convection below the neutrino sphere, and 
Appendix D introduces a simple mixing scheme by which we 
achieved to reproduce in ID simulations the main effects of 
proto-neutron star convection as observed in our 2D models. 


2. One-dimensional models 

We have chosen a total of nine progenitors from different 
groups doing stellar evolution modeling. Details of the mod¬ 
els can be found in Appendix A. These models are listed there 
in Table A. 1 and cover a zero age main sequence mass (ZAMS 
mass) range from 11 to 25 Mq. They represent various types of 
pre-collapse stellar structures. 

The one-dimensional core-collapse simulations were per¬ 
formed with our ID neutrino-hydrodynamics code Vertex with 
spectral neutrino transport, using spherical coordinates and the 
physics described in detail in Buras et al. (2006). This includes 
a state-of-the-art treatment of neutrino interactions and an ap¬ 
proximative description of general relativistic effects. The sim¬ 
ulations described in the previous and in the present paper were 
done with the equation of state (EoS) of Lattimer & Swesty 
(1991), supplemented by a general lepton-photon-baryon EoS 
(ideal gases with Coulomb corrections) that extends to den¬ 
sities below those described by the Lattimer & Swesty EoS^ 
The numerical resolution used for hydrodynamics and neutrino 
transport was also specified in Buras et al. (2006). 

The evolution can be separated into the phases of collapse, 
bounce, prompt shock propagation, neutrino burst, and accre¬ 
tion phase, which is in some cases accompanied by a transient 
shock expansion. None of our ID simulations yields a prompt 
or delayed explosion. 

The phases of core collapse, shock formation and propaga¬ 
tion, and Ve burst at shock breakout reveal only little differences 
between the progenitors because the core structure and proper¬ 
ties are very similar or become very similar during collapse (for 
details, see Appendix B and also Liebendorfer et al. 2002). At 

^ The error in the cr-particle mass fraction that was recently dis- 
coved in the EoS of Lattimer & Swesty (1991) was tested to have 
no important influence on our simulations; for details see Buras et al. 
(2006). In some models, we replaced the Lattimer & Swesty EoS be¬ 
low a density of 10^^ gcm“^, i.e. in the regime of possible cr-particle 
presence after bounce, by our low-density EoS without discovering 
any significant dynamical differences. 


the moment when the Vg burst is emitted the shock has reached 
an enclosed mass of ~ 1.0 Mq outside of which differences 
in the progenitor structure become larger. Therefore the mass 
infall rate through the shock begins to differ between the mod¬ 
els. At the same time the postshock velocities become negative 
and the shock has converted to a stalled accretion shock. Since 
the rate of mass accretion by the shock is very high shortly 
after bounce, the shock is still pushed outward for some time 
due to matter accumulating behind it. This “passive” expan¬ 
sion of the shock comes to an end when the mass accretion 
rate has dropped sufficiently and the neutrino cooling and set¬ 
tling of the accreted matter withdraw support from the shock. 
The conditions during this phase change only slowly. The sit¬ 
uation can therefore roughly be characterized by steady-state 
conditions that depend on a number of parameters governing 
the accretion phase, i.e. the mass accretion rate through the 
shock, (9^Msh, the mass Mpns and radius rpNs of the proto¬ 
neutron star (PNS)^, and the neutrino luminosity Ly (for a de¬ 
tailed discussion, see Appendix B). In recent analytic studies 
(Janka 2001; Arcones Segovia 2003) hydrostatic solutions are 
discussed which fairly well describe the nearly stationary ac¬ 
cretion of the stalled shock as seen in the numerical simula¬ 
tions at later post-bounce times, provided that the model does 
not encounter quick changes such as a sudden drop of dtM^h, 
or the onset of the explosion. 

All ID models fail to explode and the shock retreats to < 
100 km towards the end of our simulations (see top panel in 
Eig. B.4 in Appendix B). What is the effect of neutrino heating 
behind the shock? To discuss this question we consider two 
timescales, following Janka & Keil (1998), Janka et al. (2001), 
and Thompson et al. (2005). The advection timescale of matter 
falling inward from the shock to the gain radius is 


^adv(0 ~ 


rr^hit) 

^^gain(0 


1 

Mr, t) 


dr , 


( 1 ) 


where the gain radius rgam is defined as the innermost radial 
position where neutrino heating dominates neutrino cooling, 
and Vr(r, t) is the radial velocity (which is negative when matter 
is accreted by the PNS). The heating timescale is 


^heat(0 




Mt) 

where Q is the net neutrino heating rate per unit of volume and 
f) = 0, (3) 


( 2 ) 


is the so-called local specific binding energy already introduced 
in Buras et al. (2006). Here, s = eint + is the specific en¬ 
ergy, eint is the internal energy per unit of mass, and v is the 
absolute value of the fluid velocity. The gravitational potential 

^ For simplicity, we identify the PNS mass and radius with the 
mass and radius enclosed by the eleetron neutrinosphere, My^ and 
respectively. Throughout this paper, neutrinospheres are defined 
as “transport spheres”, using the opacities for momentum transfer be¬ 
tween neutrinos and stellar medium in caleulating optical depths, as 
given in Eq. (28) of Buras et al. (2006). 
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Fig. 1. The advection timescales as defined in Eq. (1) versus post- Fig. 2. The timescale ratios, Tadv/Theat, as functions of time. The lines 

bounce time. The lines are smoothed over time intervals of 5 ms. are smoothed over time intervals of 5 ms. 


(^enclosed0 is Calculated taking into account only the mass en¬ 
closed by the radius r. 

While the advection timescale represents the time matter 
spends in the gain layer (between shock and gain radius), the 
heating timescale measures the time needed for neutrino heat¬ 
ing to deposit an energy equivalent to the binding energy of 
the matter. Clearly, heating is of no importance as long as 
Tadv Theat- Thus, for obtaining a neutrino-driven explosion, 
the condition Theat ~ ^adv must hold for longer than a time inter¬ 
val Theat (see the discussion in Janka & Keil 1998, Janka et al. 
2001, Janka 2001, and Thompson et al. 2005). Note that this 
condition is not necessarily sufficient for an explosion but it 
tells one when a visible shock expansion can be expected. 

Looking at the timescales and their ratio (Figs. 1, 2), one 
recognizes that most models have long heating timescales and 
thus neutrino heating is inefficient in causing shock expansion. 
The ratio Tadv/^heat is always less than 1/2 except in the two 
low-mass models, si 1.2 and nl3, where neutrino heating is 
stronger during a short period of time. 

The ratio becomes largest during the transient shock ex¬ 
pansion associated with the times when composition interfaces 
reach the shock and the entropy makes a jump. This phe¬ 
nomenon, which is most extreme at the edges of the small iron 
cores in Model nl3 at ^pb ~ 30 ms and in Model si 1.2 around 
^pb ~ 100 ms, but in a weaker form is also present in Models 
sib, s20.0, and sl5s7b2 at around ^pb = 120 ms, 135 ms, and 
170 ms, respectively, occurs because the density and therefore 
the mass accretion rate of the infalling matter drops at the inter¬ 
faces of shells of different composition (Appendix A). The sud¬ 
den decrease of the ram pressure allows for a transient shock 
expansion until the shock finds a new equilibrium position at a 
larger radius. 


With a larger shock radius the postshock velocities are 
lower and the advection timescale increases. More efficient 
neutrino heating further strengthens the shock, the shock can 
therefore expand to even larger radii. This behaviour is ob¬ 
tained in Model nl3, and especially in Model si 1.2, where 
Tadv/^heat IS close to Unity for a period of about 20 ms. This, 
however, is shorter than the heating timescale Theat - 30 ms. 
The pronounced growth of the shock radius and of the advec¬ 
tion timescale for Model si 1.2, which results in the large local 
maxima of the corresponding curves in Figs. 1, 2, and B.4, is 
produced by a continuous strong decrease of the mass accre¬ 
tion rate, see Fig. B.4. After 113 ms post bounce this phase 
is over and the mass accretion rate continues to decline less 
rapidly. The shock then retreats quickly and finds a new quasi¬ 
stationary radius, where, however, neutrino heating becomes 
less efficient again. 

In Models sl5s7b2, s20.0, and sib the drop of dtM^^ leads 
also to shock expansion, but the effect is not strong enough to 
change Tadv significantly; also, in these models the composi¬ 
tion interface reaches the shock so late that the shock is already 
deep in the gravitational potential well and the postshock ve¬ 
locities are considerably higher than in Models nl3 and si 1.2. 
Therefore Tadv/^^heat remains well below unity. 

In summary, we find that all our ID simulations evolve 
both qualitatively and quantitatively in a similar way. In spite 
of maximum shock radii around 130-150 km the models do 
not reveal explosions. Only in the two lightest progenitor mod¬ 
els, si 1.2 and nl3, the drop in dtM^^ at the composition shell 
interfaces is sufficiently steep and large and happens suffi¬ 
ciently early to allow the shocks to reach a radius of 170 km. 
Nevertheless, the models are far from producing explosions be¬ 
cause the advection timescales remain always shorter than the 
timescales for neutrino heating, and the phases where the ra- 
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Fig. 3. Number of e-foldings that the amplitude of perturbations is es¬ 
timated to grow during advection of the flow from the shock to the 
gain radius in our ID models (cf. Eq. 4). The lines are smoothed over 
time intervals of 5 ms. For Model nl3, Wgrow starts with values around 
20 at tpb = 25 ms, and begins to decrease at tpb = 35 ms. Note that 
the entropy profile was smoothed before calculating the Brunt-Vaisala 
frequency, see the discussion in Buras et al. (2006, Section 2.5). The 
dotted line marks the instability criterion in advective fiows according 
to Foglizzo et al. (2005), and corresponds to an amplification factor of 
about 20. 


It predicts instability in static layers if Cl > 0. The Brunt- 
Vaisala frequency denotes the growth rate of fluctuations, if 
it is positive (instability), and the negative of the oscillation 
frequency of stable modes, if it is negative. In the gain layer 
between shock and gain radius, however, the gas is falling in¬ 
ward and the instability condition for Ledoux convection is not 
just given by Cl > 0 (Foglizzo et al. 2005). Here the param¬ 
eter ^grow is of crucial importance and represents the number 
of e-foldings which short-wavelength perturbations will tran¬ 
siently be amplified during their advection from the shock to 
the gain radius. Since advection has a stabilizing influence, the 
threshold for convective instability in the gain layer is found to 
be ^grow ~ 3, and the growth of modes of lowest order (i.e., 
longest wavelengths) becomes possible only when ^gi-ow ^ 5-7 
(Foglizzo et al. 2005). Figure 3 shows that only some of our 
models get close to the critical value of ^grow for which Ledoux 
convection can be expected according to a linear perturbation 
analysis. In case of Models si 1.2 and nl3, the value of ^grow 
clearly exceeds the critical threshold because in these models 
high entropy jumps at shell interfaces cause the shock to reach 
particularly large maximum radii (see above in Fig. B.4). This 
means that the postshock velocities entering the denominator in 
the integrand of Eq. (4) become smaller than in other models. It 
should be noted, however, that the discussion in Foglizzo et al. 
(2005) applies exactly only when the initial perturbations are 
very small. In case of our relatively large initial perturbations 
of order 1%, the fluctuations can grow to the non-linear regime 
— which is not accessible to the stability analysis — already 
for smaller values of ^grow 

3. Two-dimensional models 


tio of both approaches unity are much shorter than the heat¬ 
ing timescale itself. Therefore neutrino heating is not strong 
enough to drive an explosion. 

The one-dimensional models analyzed so far have one sig- 
niflcant shortcoming: They do not take into account hydrody¬ 
namic instabilities in the stellar core. Convection, especially the 
so-called hot bubble (HB) convection in the gain layer, has been 
seen to strengthen the shock in previous multi-dimensional su¬ 
pernova simulations. We can analyze our ID models for the 
existence of Ledoux-unstable regions. For this purpose we in¬ 
troduce the new variable (deflned as parameter x in Foglizzo 
et al. 2005), 

^grow(0= (4) 

JrgainCO ^r{r , t) 

where Vr is the radial component of the velocity and the Brunt- 
Vaisala frequency is deflned as 


= sgn(CL(r',r)) 




p{r',f) 


dr 


(5) 


with d0 /dr being the local gravitational acceleration. Cl is the 
Ledoux-criterion, which is given by 


Cl = 


\ds 


dv I ^p\ 

^Y^ pdr d'- 


( 6 ) 


For our two-dimensional studies with assumed azimuthal sym¬ 
metry we used the numerical code, input physics, and resolu¬ 
tion as described in Buras et al. (2006) and specifled in Table 1. 
The innermost core of 1.7 km radius was computed in spheri¬ 
cal symmetry. We have chosen three representative progenitors: 
si 1.2, which shows favorable conditions for developing strong 
Ledoux convection in the neutrino-heated postshock layer, as 
well as the two less promising progenitors, sl5s7b2 and s20.0. 
A total of seven simulations were performed in 2D, see Table 1. 
Most simulations were started from a ID model around 7 ms 
after bounce, at which time the radial velocity was randomly 
perturbed with an amplitude of +1%. For each progenitor we 
calculated a model with low angular resolution (2.7°) with a 
86.4° lateral wedge around the equatorial plane of the polar 
grid. In addition, we calculated high-resolution (1.35°) versions 
for the two lighter progenitors si 1.2 and sl5s7b2. The corre¬ 
sponding 15 Mq 2D simulation was started at the onset of core 
collapse with an initial random density perturbation (in Model 
sl5_64_p) of ±2% in order to address the question whether the 
onset of convection after bounce changes when nonradial per¬ 
turbations in the stellar core are followed through infall instead 
of being imposed shortly after bounce on a ID collapse model 
in mapping the latter onto a 2D grid. Another simulation for 
the 11.2 Mq progenitor. Model sll2_128_f, with a resolution 
of 1.41°, was performed with a full 180° grid. Finally, one sim¬ 
ulation with the 15 Mq progenitor. Model sl5_64_r, included 
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Table 1. Parameters of computed 2D models for progenitor stars with different masses. Qi is the angular velocity of the Fe-core, which is 
assumed to rotate uniformly, prior to collapse, Oq and 6i are the polar angles of the lateral grid boundaries, and No is the number of grid points 
in the lateral direction. ^20 denotes the time (relative to the bounce time) when the simulation was started or continued in 2D. 


Model * 

Progenitor 

Qi 

(rad s"^) 

[ 60 , 61 ] 

(degrees) 

No 

Resolution 

(degrees) 

Collapse 
in 2D 

Ld 

[ms] 

Perturbation 

(%) 

Boundary 

eonditions 

sll2_32 

sll.2 

- 

[46.8,133.2] 

32 

2.70 

- 

6.4 

V,+l 

periodie 

sll2_64 

sll.2 

- 

[46.8,133.2] 

64 

1.35 

- 

7.3 

V,+l 

periodic 

sll2_128_f 

sll.2 

- 

[0,180] 

128 

1.41 

- 

6.4 

V,+1 

reflecting 

sl5_32 

sl5s7b2 

- 

[46.8,133.2] 

32 

2.70 

- 

6.5 

v,+l 

periodic 

sl5_64_p 

sl5s7b2 

- 

[46.8,133.2] 

64 

1.35 


-175 

P, +2 

periodic 

sl5_64_r 

sl5s7b2 

0.5 

[0,90] 

64 

1.41 

- 1 - 

-175 

P,+1 

reflecting 

s20_32 

s20.0 

- 

[46.8,133.2] 

32 

2.70 

- 

6.5 

V,+1 

periodic 

sl5_mix ^ 

sl5s7b2 

- 

- 

1 

- 

- 


- 

- 


* Models sll2_32, sll2_64, sl5_32, sl5_64_p, sl5_64_r, and s20_32 were discussed in overview in Buras et al. (2003). 

^ This is Model sl5Gio_32.b from Buras et al. (2006). Different from the other models presented in this paper. Model sl5_32 was computed 
with an older implementation of the gravitational potential, which contained a slightly different treatment of the relativistic corrections in the 
potential (see also the comment in Table B.l). This affects only the neutrino luminosities and average neutrino energies by a few percent. We 
have corrected this difference in our plots for the luminosities and mean energies. The correction was calculated as AL = LiD,new - TiD,oid. where 
^iD,new and LiD,oid arc the luminosities of the ID simulations with the “new” and “old” versions of the gravitational potential, respectively. Then 
the corrected luminosity is given by L 2 D,corr = T 2 D,oid + AL, where L 2 D,oid is the luminosity of the 2D model which was simulated with the old 
version of the gravitational potential. For the average neutrino energies the procedure is analogous. 

^ This model was calculated in ID using the mixing algorithm described in Appendix D for treating effects from PNS convection. 


rotation. The angular frequency at the onset of core collapse 
was assumed to be Qi = 0.5 rad s“^ and constant in the Fe and 
Si core, and decreasing (spherically symmetrically) like 
outside of 1750 km (1.43 Mq). This choice of the rotation rate 
and rotational profile was motivated by results for pre-collapse 
stellar cores of stars whose evolution is followed including the 
angular momentum transport by magnetic fields (Heger et al. 
2005). Our choice of Qi in the iron core is roughly ten times 
lower than the core rotation of non-magnetic stars (Woosley 
et al. 2002) and about ten times larger than Heger et al.’s rota¬ 
tion rates of magnetized stars. It intends to maximize the effects 
of rotation during core collapse under the constraints that (a) 
the initial star can well be considered as spherically symmetric, 
and that (b) the assumed rotation, which is superimposed on a 
ID stellar model, does not imply significant deviations from the 
hydrostatic equilibrium and gravitationally bound state of the 
progenitor model. In order to fulfill both requirements we limit 
the rate of rotation such that the influence of centrifugal forces 
is very small prior to collapse: F’cent/^grav <1% everywhere 
for the initial ratio of the centrifugal force to the gravitational 
force. The density of the stellar model was perturbed by +1% 
and the simulation was carried out with a lateral resolution of 
1.41° in a wedge from the polar axis to the equator (i.e. besides 
axial symmetry also equatorial symmetry was assumed). 

In Sect. 3.1 we shall focus on convection below the neu- 
trinosphere of the proto-neutron star (“PNS convection”). We 
will try to classify this convection and will describe its effects 
on the evolution of the supernova and on the neutrino emission. 
To this end we will develop a mixing algorithm as described in 
Appendix D and will present a ID simulation performed with 
it. We shall finish by discussing resolution and perturbation 
effects and differences between the simulations with different 
progenitors. In Sect. 3.2 a description of convective overturn in 
the neutrino-heating layer behind the shock (“hot bubble (HB) 


convection”) will follow. Again, we will first discuss differ¬ 
ences compared to ID models. Here, the resolution and size 
of the angular wedge play a much more important role and we 
need to elaborate on both aspects. Also the sensitivity of the on¬ 
set of HB convection to the size of the seed perturbations will 
be discussed, and our results for the evolution of pre-collapse 
perturbations during the infall phase in Model sl5_64_p can be 
found in Appendix E. Since differences of this model compared 
to Model sl5_32 were minor, the results of this simulation will 
be included in the plots following below mostly without any 
special discussion. In Sect. 3.3 we shall describe our calcula¬ 
tion with the full 180° grid, which allows low (I = 1,2) modes 
to develop with important consequences. Section 3.4 will con¬ 
tain our results for the rotating 15 Mq model. We shall discuss 
the influence of rotation on the PNS, its neutrino emission, and 
convection according to the Solberg-Hpiland-criterion (gener¬ 
alized to include the effects of Fe-gradients on convection). 
The role of centrifugal forces for the HB convection and shock 
propagation will also be addressed. 


3.1. Proto-neutron star convection 

Many phenomena inside the PNS, e.g. the conditions for con¬ 
vection in its interior, can mostly be discussed without taking 
into account the convective activity in the HB layer (but not in¬ 
versely!). For several of the discussed models (sl5_32, s20_32, 
sl5_64_p) this is true, because they develop only weak HB con¬ 
vection with little impact on the shock propagation and on the 
PNS. Therefore nearly steady-state conditions prevail around 
the PNS, the mass in the gain layer changes only very slowly, 
and the mass accretion rate onto the PNS is approximately the 
same as the mass accretion rate through the shock. Moreover, 
the accretion flow onto the PNS is nearly laminar and isotropic. 
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Stability criterion In order to analyze the hydrostatic neutron 
star for convective instability we consider the Brunt-Vaisala 
frequency, Eq. (5). However, since neutrinos are strongly cou¬ 
pled with the dense plasma and close to equilibrium with 
the matter in the PNS, we generalize the Ledoux-criterion 
to a “Quasi-Ledoux criterion”, in which the effects of neu¬ 
trino transport are approximately accounted for (see Wilson & 
Mayle 1993, Buras et al. 2006), 


Cql = 


dp \ d (vtot) 

dp 




dr 

^d(F,ep) 


'®p ’{sMp) 


dr 


■^difF 


dFi, 


iep 


dr 


> 0 , 


(7) 


for instability, where Vtot = ^ is the entropy including the 
neutrino entropy, Fiep = Fe + is the electron lepton number, 
and y^difF ~ 1 is a parameter which expresses the efficiency of 
lepton number exchange between buoyant fluid elements and 
their surroundings via neutrino diffusion. For our evaluation we 
choosey^difF = 1, following Keil (1997). The brackets () denote 
angular averages, whereas the last term,ySdifFdFiep/dr, involves 
local evaluation because it describes the change of lepton num¬ 
ber along the paths of buoyant fluid elements by the neutrino 
exchange between such fluid elements and their surroundings. 
Note that due to the nearly perfectly central gravitational ac¬ 
celeration in case of non-rotating neutron stars, only the radial 
components of the gradients of Vtot and Fiep occurring in the 
general vector form of the Ledoux criterion (cf. Eq. 9) are rel¬ 
evant (buoyancy acts only perpendicular to equipotential sur¬ 
faces). Moreover, we simplified the analysis by performing it 
on the laterally averaged stellar background. The determination 
of stable and unstable layers on the basis of laterally averaged 
stellar quantities leads to essentially the same conclusions as 
the lateral average of the Ledoux criterion after evaluation with 
the 2D data. 

Recently, Bruenn et al. (2004) presented an elaborate dis¬ 
cussion of hydrodynamic instabilities in the PNS including the 
effects of neutrino diffusion (an extension of a previous anal¬ 
ysis by Bruenn & Dineva 1996). They argue that local pertur¬ 
bations in the lepton number will be reflected in the neutrino 
phase space and thus cause a net neutrino diffusion which tries 
to wash out the perturbations, an effect which can be accounted 
for by a “response function”. Since neutrinos also carry en¬ 
tropy, the neutrino diffusion that smoothes the lepton number 
perturbations will create entropy perturbations. This effect is 
characterized by a “cross response function”. Of course, en¬ 
tropy perturbations will in an analogous way induce an equi¬ 
librating net neutrino diffusion which at the same time trans¬ 
ports lepton number between the fluid elements and their sur¬ 
roundings, corresponding to another “response function” and 
a “cross response function”. Bruenn et al. (2004) found in a 
numerical analysis that perturbation-induced neutrino lepton 
number transport by diffusion is considerably more rapid than 
thermal transport, and that the transport of lepton number reacts 
faster to entropy perturbations than to lepton number perturba¬ 
tions. For such a situation convective instability should set in 
for most stellar conditions, even when the fluid is Ledoux sta¬ 
ble. In particular, Bruenn et al. (2004) describe two kinds of 



r [km] 



e [MeV] 

Fig. 4. Top: Transport mean free paths of muon and tau neutrinos and 
antineutrinos for different energies as funetions of radius at 243 ms 
after bounce, about 200 ms after the onset of PNS convection in Model 
sl5_32. Also shown is the local density scale height (bold solid line). 
The vertical line marks the radius of the outer edge of the convective 
layer in the PNS (cf. Fig. 5). Bottom: Mean free paths of Ve, Ve^ and 
as functions of neutrino energy at a radius r = 25 km (marked by the 
vertical line in the upper panel). The arrows indicate the mean energies 
of the Ve flux (solid), Ve flux (dashed), and flux (dotted, coinciding 
with the dashed arrow) at this position in the star. 


instabilities in the presence of neutrino diffusion: One insta¬ 
bility occurs when the entropy difference between a displaced 
fluid element and its surroundings in a background with an en¬ 
tropy gradient results in a lepton fraction difference, which pro¬ 
vides a driving force such that the induced perturbation grows. 
The buoyant rise of a perturbed fluid element, together with 
neutrino diffusion, will thus further increase the difference in 
entropy between the fluid elements and their surroundings and 
will create lepton number fluctuations from entropy perturba¬ 
tions, which continue to drive buoyant motion. A second in¬ 
stability exists where the neutrino diffusion creates an “over¬ 
stable” situation, i.e. where the effect of neutrino diffusion will 
drive a perturbed fluid element back to its original position, but 
to such an extent that the fluid element overshoots and thus 
oscillates around its original position with increasing ampli¬ 
tude. Bruenn et al. (2004) call these two doubly-diffusive in¬ 
stabilities “lepto-entropy Anger” (LEE) convection and “lepto- 
entropy semiconvection” (LESC), respectively. They also dis¬ 
tinguish Ledoux convection. However, Ledoux and LEE con¬ 
vection are closely related (LEE convection is an extension of 
Ledoux convection). 

Bruenn et al.’s analysis of stellar profiles shows that 
Ledoux/LEF convection should appear in an extended region 
of the PNS from around 15 km to the neutrinospheres, whereas 
LESC should be visible deeper in the core. We suppose that our 
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Fig. 5. Snapshots of PNS convection in Model sl5_32 at 48 ms (a) and 243 ms (b) after bounce. The upper left quadrants of each plot depict 
color-coded the absolute value of the matter velocity, the other three quadrants show for Ve, Ve^ and Vx (clockwise) the ratio of the local neutrino 
flux to the angle-averaged flux as measured by an observer at inflnity. The diagonal black lines mark the equatorial plane of the computational 
grid and the thick circular black lines denote the neutrinospheres (which have a radius larger than 60 km in flgure a). The flgures (c) and (d) 
show at the same times the relative lateral variations of lepton number and entropy (including neutrino entropy), i.e. SX = (X - (X)^)/ (X)^ for 
a quantity X. 


stellar profiles might yield qualitatively similar results if one 
applied their analysis. Our two-dimensional simulations, how¬ 
ever, do not confirm their predictions in detail. For example, 
Model sl5_32 shows convective instability between 17 km and 
30 km at ^pb = 48 ms, see Fig. 5a. We interpret this as Ledoux 
convection and determine growth rates of the Ledoux instabil¬ 
ity (cf. Figs. 6 and 9) that are typically more than a factor of 10 
higher than those of LEF found by Bruenn et al. (2004). Our 
hydrodynamic simulations reveal a behavior in agreement with 
our larger rate estimates. Moreover, the region showing con¬ 
vective activity does not extend to the neutrinospheres, which 
are beyond 60 km at that time (Fig. 11), and also at later times 
the upper boundary of the convective layer stays well below the 
neutrinospheres (Fig. 5b). 

We suspect that the discussion by Bruenn et al. (2004) is 
not fully applicable up to the neutrinospheres. It is based on 
the assumptions that the mean free path (mfp) of neutrinos is 
much smaller than the size of the displaced fiuid elements and 
shorter than local gradients in the background medium so that 
neutrinos diffuse and stay close to local equilibrium. In con¬ 
trast, Fig. 4 reveals that the neutrino mfp becomes large near 
the upper edge of the convective zone, in particular for muon 
and tau neutrinos. As a consequence, the neutrinos start stream¬ 
ing increasingly rapidly and the neutrino densities in this region 
(first for and Vj, then for Vq and finally for Ve) begin to de¬ 


viate from local chemical equilibrium. The net rate of neutrino 
losses therefore grows strongly, leading to a quick rise of the 
neutrino luminosities. Despite their growing mfp, the produc¬ 
tion of neutrinos in the stellar medium is sufficiently strong to 
allow for efficient energy and lepton number drain. We there¬ 
fore suspect that in the semitransparent layer below the neutri- 
nosphere the accelerating radial transport causes enhanced neu¬ 
trino release such that entropy and lepton number fluctuations 
are wiped out by neutrinos streaming off. This happens before 
local exchange between fiuid elements and their surroundings 
can efficiently take place as assumed by Bruenn et al. (2004). 
LEF instabilities thus have no time to grow. 

LESC as predicted by Bruenn et al. (2004) to occur far in¬ 
side the PNS is not visible in our models. However, these au¬ 
thors mention that the existence of LESC is very sensitive to 
the exact values of the response functions, which in turn de¬ 
pend on the details of the neutrino interactions. We suspect that 
the different description of neutrino-matter interactions in our 
simulations might prevent LESC. In any case, convection very 
far inside the PNS should have less influence on the shock dy¬ 
namics than the PNS convection which we see in our models. 
Therefore we will ignore the possibility of LESC in the follow¬ 
ing discussion and consider its effects on the neutrino emission 
and the explosion mechanism as negligible. 
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Fig. 6. Brunt-Vaisala frequency (Eq. 5), using Eq. (7) as stability cri¬ 
terion, evaluated for different ID models at 20, 30, and 40 ms after 
bounce. 


Phenomenology In Figs. 6-8, we see that a convectively un¬ 
stable layer begins to develop in the PNS about 30 ms after 
bounce. This happens in a region which initially is stable due 
to a positive entropy gradient, see Fig. 7. As neutrino diffu¬ 
sion carries entropy away more efficiently at larger radii where 
the optical depth is lower, the entropy profile flattens, and fi¬ 
nally the entropy gradient turns negative and the PNS becomes 
Ledoux unstable. Ten ms later the perturbations have grown 
(with a growth rate of about ojbw ~ 0.5-1.0ms“\ see Figs. 6 
and 9) to become non-linear and strong PNS convection has set 
in. This can be seen in the large lateral velocities in Fig. 8 (dark 
shaded). Similar to Keil (1997), we find that the structure of the 
convective cells is initially that of rolls with angular sizes be¬ 
tween 20° and 30° and radial extension between 10 and 15 km 
(Fig. 5a). The rolls are stable for about 5 ms, then decay and 
form again at different locations in the same layer. 200 ms later, 
the contraction of the PNS has reduced the radial size of the 
convective cells to 10 km. The overturn velocities are around 
3 X 10^ cms“^ with peaks of up to 5 x 10^ cms“^ At later 
times the velocities decrease to values around 2 x 10^ cms“^ 
The region with large lateral velocities (>700 km s“\ dark 
shaded in Fig. 8) is wider than the layer with instability ac¬ 
cording to Eq. (7). This has two reasons: First, rising or sinking 
fluid elements can over- and undershoot into the adjacent sta¬ 
ble layers. Second, the periodic boundary conditions applied in 
some of our simulations allow rings of uniform, lateral velocity 
to occur. This artificial phenomenon is associated with matter 
that settles from the neutrino-heated convective postshock layer 
onto the PNS and has obtained large lateral velocity compo¬ 
nents by participating in the overturn motion in the gain layer. 
The rings show up in Fig. 8 as a layer with high lateral veloci¬ 
ties that moves from the “hot-bubble” (HB) convective zone to 
the convection zone in the PNS between 100 and 150 ms post¬ 
bounce. Looking at the velocity distribution in 2D snapshots 



Fig. 7. Lepton number and total entropy versus enclosed mass in the 
PNS for the ID model sl5s7b2 for different post-bounce times before 
the onset of PNS convection in the corresponding 2D simulation. 
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Fig. 8. Convective region in Model sl5_32. The dark-shaded regions 
have lateral velocities above 7 x 10^ cms"\ the light-shaded regions 
fulfill the Quasi-Ledoux criterion, Eq. (7). The solid lines mark the 
radii up to 30 km in 5 km steps, and 50 km; the dashed lines indicate 
density contours for 10^"^, 10^^, 10^^, and 10^^ gcm“^. 


of the simulation. Fig. 5b, the rings can be identified around 
r=30 km. 
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Fig. 9. Brunt-Vaisala frequency in the ID simulation of the stellar pro¬ 
genitor sl5s7b2 (thin) and the corresponding 2D simulation, Model 
sl5_32 (thick), at 30 ms (dashed), 62 ms (solid), and 200 ms (dash- 
dotted) post-bounce. At 30 ms, the thin and thick dashed lines coin¬ 
cide. All lines are truncated at 50 km. Positive values indicate insta¬ 
bility. For the 2D model the evaluation was performed with laterally 
averaged quantities. 

Throughout our simulations, PNS convection occurs exte¬ 
rior to an enclosed mass of 0.6 Mq^, see Fig. 8. This inner 
boundary changes only little during the simulations, whereas 
the outer boundary of the convective layer moves outward in 
mass as time goes on, following the ongoing accretion of mat¬ 
ter on the PNS. Keil et al. (1996) found in their models that 
PNS convection develops in an initially narrow layer, but the in¬ 
ner edge of this layer moves continuously deeper into the PNS. 
Their models, however, were evolved until 1.3 s after bounce, 
and the inward motion of the lower boundary of the convective 
zone may take place at times not covered by our present cal¬ 
culations. Moreover, the velocity at which the convective layer 
digs deeper into the PNS must be expected to depend on the 
treatment of the neutrino transport, which was described very 
approximatively by grey, radial, flux-limited equilibrium diffu¬ 
sion by Keil et al. (1996) and Keil (1997), who also made quite 
radical approximations for the neutrino-matter interactions (see 
Appendix in Keil & Janka 1995). 

We note here that a closer discussion of the first 200- 
300 ms of post-bounce evolution, comparing these older diffu¬ 
sion simulations with our present models, is not useful to con¬ 
vey our understanding particularly of the possible deficiencies 
of the simpler diffusion treatment. Such a comparison is ham¬ 
pered by other very important differences between the simula¬ 
tions. Firstly, the post-collapse models from which the calcu¬ 
lations of Keil et al. (1996) and Keil (1997) were started, were 
provided by S. Bruenn and thus were computed with a neu- 

^ Note that in 2D simulations the “enclosed mass” is not a 
Lagrangian coordinate but corresponds to the sphere which contains a 
certain amount of mass at a given time. 


trino transport treatment that was not compatible with the one 
used for the subsequent 2D evolution. This must be suspected 
to have caused transient effects during an early phase of the 
post-bounce evolution of undetermined length. Such transients 
are not present in our current, consistent models. Secondly, the 
simulations by Keil et al. (1996) and Keil (1997) considered 
a neutron star of l.l-1.2Mo and thus much smaller in mass 
than the compact remnants of the models in this paper, which 
grow in mass particularly rapidly during the first hundreds of 
milliseconds after bounce when the mass accretion rates are 
still high (see Fig. B.4). And thirdly, neglecting the accretion 
layer outside of the neutron star and placing the outer grid 
boundary near the neutrino sphere by Keil et al. (1996) and Keil 
(1997) may also have caused differences of their results from 
the present, full simulations of the collapsing supernova core. 
We therefore refrain from attempting a closer comparison of 
these old and our new models beyond the level of qualitative 
statements about the existence, long-lasting presence, and ba¬ 
sic structural features of a convective layer inside the nascent 
neutron star. 

The convective flow transports energy and lepton number 
from deeper layers of the PNS closer to the neutrinospheres, 
thereby flattening the entropy gradient and the lepton number 
profile, see Fig. 10. However, at the end of the calculation the 
initially quite steep negative Fiep gradient has not disappeared 
completely. This suggests that the convective mass motions can 
not efficiently transport lepton number over large distances. 
The reason for this is the fact that buoyant mass elements ex¬ 
change leptons easily with their local surroundings via neutrino 
diffusion. Therefore rising bubbles with initially large Fiep end 
up with low Fiep when they reach the outer stable layers of the 
PNS. For this to happen the timescale of Fiep-equilibration be¬ 
tween buoyant bubbles and their surroundings via neutrino dif¬ 
fusion must be of the same order as the rise timescale of the 
bubbles. This situation corresponds to a value of y^diff in Eq. (7) 
close to unity. If the rise time were much longer, the ID and 2D 
profiles of Fiep in the convective region would not differ, while 
if it were much shorter, the Fiep profile would be flat like the 
entropy profile. 

A second effect contributes in establishing the observed 
conditions in the convective PNS layer. In contrast to the energy 
loss from the stellar medium, which is small in the considered 
period of time compared to the huge heat capacity of the accre¬ 
tion layer, lepton number is radiated away very efficiently, at 
least initially. At early times, when Fiep is still large in the con¬ 
vective PNS layer, the number flux of Vq is much larger than 
that of Ve and the lepton number loss proceeds very rapidly. At 
later times, however, it slows down because Fiep has decreased 
and the release of Ve and Vg has become more similar. 

Effects on radial structure and neutrino emission of the PNS 
The redistribution of lepton number and energy gradually af¬ 
fects the structure of the PNS. The “drain” region, where Fiep 
and internal energy eint are extracted, has densities close to the 
nuclear density (0.3-2 x 10^^ g cm“^), whereas the “dump” re¬ 
gion, where lepton number and energy are deposited, has lower 
densities (lO^^-lO^^ gcm“^), see Fig. 8. A detailed analysis 
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Fig. 10. Lepton number, total entropy, specific internal energy, and 
density profiles versus enclosed mass in the PNS for a sample of ID 
(dotted) and 2D simulations of different progenitor stars, 200 ms after 
bounce (i.e., approximately 160 ms after the onset of PNS convection). 
For the 2D models angle-averaged quantities are plotted. 

presented in Appendix C reveals that the transport of lepton 
number from the drain region to the dump region will effec¬ 
tively lead to expansion. Also the transfer of energy from the 
drain region to the dump layer leads to an expansion of the 
PNS. Consequently the neutrino sphere radii in 2D simulations 
are larger than in the corresponding ID models, see Fig. 11. 

The changes in the radial profiles of energy and lepton 
number and the expansion of the PNS relative to ID simula¬ 
tions have interesting consequences for the neutrino emission. 
Convection in the PNS alters the emission of neutrinos in and 
above the dump region and therefore the PNS loses energy and 
lepton number at a different rate than without convection. In 
this case the linear analysis (Appendix C) shows that an in¬ 
creased loss of lepton number in and above the dump region 



tpb [ms] 

Fig. 11. Radius of the electron neutrinosphere (as a measure of the 
PNS radius) for the 2D models and their corresponding ID models. 
For Model sl5_64_r we show the equatorial (upper) and polar (lower) 
neutrinospheric radii. For other 2D models, angle-averaged quantities 
are plotted. 


leads to a more extended PNS, while an increased loss of en¬ 
ergy leads to a more compact PNS (Fig. C.l). As will be dis¬ 
cussed in detail below, PNS convection causes enhanced lepton 
number release and initially reduced energy loss via neutrinos, 
both supporting the “expansion” of the PNS relative to the ID 
simulations. At times later than 150ms post bounce, the energy 
loss in neutrinos is also enhanced relative to the ID models, 
without however being able to compensate the inflative effects 
on the PNS structure caused by the convective redistribution of 
energy and lepton number and by the enhanced lepton losses 
form the neutrinospheric layer due to neutrino emission. 

Immediately after the onset of PNS convection the increase 
of Fiep near the upper boundary of the convective region leads 
to a higher electron chemical potential /ie and therefore a higher 
electron neutrino chemical potential= /ie +yUp -/in. Thus the 
Ve abundance (given by equilibrium conditions) decreases sig¬ 
nificantly in this layer. Although this happens in a region where 
the Ve luminosity has reached only 5-10% of its final value (see 
Fig. 5a) the decrease is sufficiently large to lower the total Vq lu¬ 
minosity by several percent. When the PNS convection has de¬ 
veloped to full power, Ly^ can be lowered by up to 10% in com¬ 
parison to the ID models, see Fig. 12. The decrease of the Vq 
equilibrium abundance also affects the heavy-lepton (“Vx” for 
Vyu, Vr, or Vj) luminosity via the process VeVe ^ VxVx, which 
is the dominant Vx production process in the region where the 
reduction of the Ve abundance happens. However, Ly^ in that 
region has already attained 70% of its final value, and there¬ 
fore the reduced Vx production rate leads to a decrease of the 
radiated Vx luminosity by at most 5%, typically less (Fig. 12). 

After 100 ms post-bounce these moderate effects are over¬ 
ridden by the structural changes of the PNS, which lead to 
larger neutrinospheric radii than in ID models, as well as the 

















12 R. Buras et al.: Two-dimensional hydrodynamic core-collapse supernova simulations with spectral neutrino transport 



0 50 100 150 200 250 

tpb [ms] 


Fig. 12. Luminosities of electron neutrinos, Ve, (top) electron an¬ 
tineutrinos, Ve. (middle), and one kind of heavy-lepton neutrinos, 
(denoted by Vx; bottom) for the 2D models and for their 
corresponding ID counterparts, evaluated at a radius of 400 km for 
an observer at rest. Note that the luminosities of Model sl5_32 were 
corrected for the differences arising from the use of a slightly different 
effective relativistic potential as described in the context of Table 1. 

ongoing convective transport of energy into the region below 
the Vx-sphere. Moreover, the convective layer now extends to 
lower densities so that the region affected by PNS convection 
contributes now 30%, 30%, and 90% to the radiated Vg, Ve. and 
Vx luminosities, respectively (see Fig. 5b), although this region 
has a transport optical depth (for a definition see Buras et al. 
2006, Eq. 28) in excess of 40, 14, and 10 for Vg, Vg, and Vx, re¬ 
spectively, and therefore is still well inside the neutrino sphere. 

Larger neutrino spheric radii without the described convec¬ 
tive inflow of energy would lead to lower luminosities as a con¬ 
sequence of an associated decrease of the neutrino spheric tem¬ 
perature Ty. This, for example, is seen in simulations with dif¬ 


ferent EoSs where the PNS radius depends on the high-density 
EoS properties. A larger PNS radius correlates with lower lu¬ 
minosities (Janka et al. 2005). In contrast, in our 2D models 
the luminosities increase. We indeed find lower Ty, which re¬ 
sult in lower mean neutrino energies {Sy) (defined by the ratio 
of energy to number flux), see Pig. 13. The difference can be up 
to 10% for all neutrino kinds after 200 ms of PNS convection. 
Because of the energy transport to the neutrino spheres by con¬ 
vection, however, this reduction in Ty is much weaker than it 
would be in an adiabatically expanding layer. Apparently, the 
larger neutrinospheric radii and only slightly lower tempera¬ 
tures lead to a net increase of the luminosities relative to the 
ID results, see Pig. 12. The effect is strongest for Vx, which 
decouple energetically from the medium already near the up¬ 
per boundary of the convective layer (Pig. 5b); after 200 ms of 
convection, Ly^ is almost 20% higher than in the ID models. 
Ly^ increases only by a few percent, while Ly^ is almost identi¬ 
cal in ID and 2D models, which means that the higher electron 
chemical potentials and the effects associated with the convec¬ 
tive energy transport and structural changes of the PNS nearly 
compensate each other. 

After the onset of PNS convection the 2D models delep- 
tonize faster than their ID counterparts (Pig. 14). The lepton 
number loss is enhanced after 250 ms of post-bounce evolution 
by typically 8-10% (compare Pig. 14 with Pig. B.5). The evolu¬ 
tion of the total energy loss is more complex and is smaller than 
in the ID simulations for the hrst ~ 100-140 ms of reduced en¬ 
ergy emission. Only afterwards the losses become stronger in 
the 2D cases. However, even 250 ms after bounce, the total en¬ 
ergy loss is only 2-4% higher with PNS convection. 

Note that the neutrinospheric luminosities in 2D models 
relative to ID models are increased by 25%, 15%, and 25% for 
Ve, Ve, and Vx after 200 ms. The effect of the PNS convection 
therefore manifests itself more strongly at the neutrinospheres 
than farther outside, because the luminosity differences get re¬ 
duced due to structural differences above the neutrinospheres 
and thus resulting differences in the neutrino absorption and 
emission outside of the neutrinospheres. 

We have performed several 2D simulations with varied res¬ 
olution (2.7° and 1.35° or 1.41°) and different choices of the 
lateral grid (either a wedge with 86.4° around the equator or 
a full 180° grid). The results show that these differences affect 
the behaviour of the PNS convection only little. The flattening 
of the Fiep prohle appears to be slightly more effective at higher 
resolution, see Pig. 10, but the difference is too small to be dis¬ 
tinguishable in other observables. Thus, a resolution of 2.7° 
seems sufficient to simulate PNS convection, mainly because 
the convective cells have a size of 20°-30° and are therefore 
much larger than the size of the angular grid zones. Similarly, 
the convective cells are sufficiently small to ht several times 
into the angular wedge in each of our simulations. Boundary 
conditions are therefore found to have no visible influence, and 
our simulation with a full 180° grid revealed no important dif¬ 
ferences of PNS convection and of its consequences compared 
to the models with a 90° angular wedge. 

The PNS convection also evolves in a similar way in our 
simulations with different progenitors. The relative differences 
compared to the ID models are quite similar for ah progenitors. 
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Fig. 13. Average energies of the radiated neutrinos (defined by the ra¬ 
tio of energy to number flux) for the 2D models and for the corre¬ 
sponding ID models, evaluated at a radius of 400 km for an observer at 
rest. The lines are smoothed over time intervals of 5 ms. Note that the 
average neutrino energies of Model sl5_32 were corrected for the dif¬ 
ferences arising from the slightly different effective relativistic gravi¬ 
tational potential as described in the context of Table 1. 

Moderate differences, e.g. in the total energy loss, f - 
f L^^dt, between Models sll_32, sl5_32, and s20_32, see 
Fig. 14, appear small considering the large differences in the 
mass accretion rates and post-bounce evolution of these stars. 
The evolution of the energy loss is clearly different and lower 
only in the rotating Model sl5_64_r, where a significant amount 
of energy is stored in rotation instead of being radiated by neu¬ 
trinos. Also different seed perturbations have little effect on the 
evolution of convection in the PNS. Large seeds as, e.g., in 
Model sl5_64_p, lead to PNS convection which starts typically 
10-20 ms earlier. This short time difference of PNS convection, 
however, has hardly any influence on the supernova evolution. 
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Fig. 14. Differences between the total lepton number (top) and energy 
losses (bottom) of the 2D models and their corresponding ID models 
as functions of post-bounce time. Here, X is the total lepton number 
flux, and L is the total neutrino luminosity. 


Effects on the gain layer The structural differences between 
ID and 2D models, which are induced by PNS convection, also 
influence the layer between the neutrinospheres and the shock. 
We can discuss these influences on the basis of Models sl5_32 
and s20_32, in which HB convection remains rather weak so 
that the differences from ID result mainly from effects associ¬ 
ated with PNS convection. In order to understand how much 
PNS convection contributes to the differences between ID and 
2D models we have also performed a ID simulation, Model 
sl5_mix, for the progenitor sl5s7b2, in which PNS convec¬ 
tion was treated with the simple mixing scheme described in 
Appendix D. In Fig. 10 we see that this approximative treat¬ 
ment reproduces the transport of energy and lepton number 
by PNS convection rather well. Therefore Model sl5_mix can 
be understood as a simulation with PNS convection but with¬ 
out HB convection. Comparing Models sl5s7b2, sl5_32, and 
sl5_mix, we conclude that model properties which basically 
are independent of the evolution of the gain layer, such as 
neutrinosphere radii and neutrino luminosities, are sensitive to 
PNS convection but not to HB convection (Figs. 11, 12). The 
differences between the ID and 2D models in the gain region, 
which are visible in the mass of the gain layer, Mgi, the ad- 
vection timescale, and the shock radius, r^h (Figs. 17, 15), 
are initially mainly caused by the HB convection, despite the 
fact that it is weak. However, the changes due to PNS convec¬ 
tion become gradually more and more important, and make the 
dominant influence at the end of the simulations. In the fol- 
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lowing we will only discuss the effects which result from PNS 
convection. Model sl5_mix will therefore also be considered 
as a “2D model”. 

Note that in two-dimensional models with downflows 
and buoyant rising bubbles in the convective, neutrino-heated 
layer, Eq. (1) is not applicable as definition of the advection 
timescale, and we therefore use the following definition of an 
“effective advection timescale”: 

^ = t2(Mi) - ti(Mi ), ( 8 ) 

where t 2 is defined by the condition M(r = rgam, t = t 2 ) = Mi 
and ti by the condition M(r = = h) = Mi. This is the 

time difference between the moment when the shock encloses 
a mass Mi and the time when this same mass M/ is enclosed by 
the gain radius. The expression is evaluated for different val¬ 
ues of Mi to monitor the time evolution of In cases where 
this definition of the advection timescale is used, we also ap¬ 
ply it to the corresponding ID models for comparison. In the 
multi-dimensional situation is a measure of the average 
storage period of gas in the gain layer, providing a rough es¬ 
timate for the time interval which accreted matter is exposed 
to neutrino heating, without however considering in detail the 
complex fluid flow in downdrafts and rising plumes. 

As explained above, PNS convection leads to larger radii of 
the neutrinospheres. The mass enclosed by the neutrinospheres, 
however, is smaller because of a less compact PNS. The struc¬ 
ture of the PNS is changed such that this also applies for the 
gain radius and its enclosed mass, as can be seen for Models 
sl5s7b2, sl5_mix, and sl5_32 at late times in Fig. 16. However, 
the mass which is enclosed by the standing shock is identical 
in ID and 2D models because it depends on the total mass ac¬ 
creted by the shock and therefore on the infall region and pro¬ 
genitor structure, which do not differ between the ID and 2D 
models (the shock radius has almost no influence on the en¬ 
closed mass because of the relatively low densities behind the 
shock). For these reasons the mass in the gain layer must be 
expected to be larger in 2D. The top panel of Fig. 17 confirms 
this conclusion. We And that the mass in the gain layer, Mgi, 
can be more than a factor of two larger than in ID. Inevitably, 
with larger gain radius and larger mass in the gain layer, the 
shock radius is also larger than in ID, see Fig. 15. 

Since the density at a given radius between the neutri¬ 
nospheres and the shock is higher in a model with PNS con¬ 
vection (this is a structural consequence of the larger Vg-sphere 
radius, which is located at p ^ 10^^ gcm“^), and the mass ac¬ 
cretion rate M{r) is approximately constant for different radii 
above the gain radius (a feature of nearly stationary conditions), 
and M is equal in ID and 2D (a consequence of the conserva¬ 
tion of the mass flow through the shock), the postshock veloci¬ 
ties are (on average) smaller in 2D. This implies a larger effec¬ 
tive advection timescale, see Fig. 17, second panel. The values 
of turned out to be larger by up to a factor of three in the 
2D simulations. 

Astonishingly, we And that the total heating rate of the gain 
layer, dtE^\, is almost identical in ID and 2D (Fig. 17, third 
panel). Model sl5_mix, however, reveals that this is a coinci¬ 
dence. A comparison of this model with Model sl5s7b2 shows 
that PNS convection actually reduces dtE^\ slightly. We will 
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Fig. 15. Maximum and minimum shock radii as functions of time for 
the 2D models of Table 1, compared to the shock radii of the corre¬ 
sponding ID models (dotted lines). 

not try to comprehend this behaviour analytically here, we only 
mention that the lower Vq and Vq luminosities and the lower neu¬ 
trino energies (which the absorption rates depend on) lead to 
weaker heating. Also the fact that the gain radius is larger may 
contribute to this decrease. We believe that the larger mass in 
the gain layer is of minor importance for StEgi, because this ad¬ 
ditional mass is located at larger radii where the heating rate is 
small. 

Since Mgi increases as a consequence of PNS convection, 
also Theat Mg\l6tEg\ increases. The ratio Tadv/^heat, however, 
changes much less (Fig. 17, fourth panel). This can be under¬ 
stood by the fact that both Mgi and Tadv increase due to PNS 
convection, as explained above. These changes partly compen¬ 
sate each other in the evaluation of the timescale ratio. 

In summary, PNS convection has two important conse¬ 
quences: First, the emitted neutrinos have lower mean ener- 
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Fig. 16. Minimum and maximum gain radii as functions of post¬ 
bounce time for our 15Mq 2D simulations compared to the gain radius 
of the corresponding ID models. 

gies than in ID models (up to 10% lower after 200 ms of PNS 
convection). The neutrino luminosities initially decrease due 
to the onset of PNS convection (Ly^ by about 10 % and Ly^ 
by about 5 %) and increase at later times {t ^ 150 ms post¬ 
bounce), to reach several percent higher values for Ve and even 
~ 20% higher values for Vx after 200 ms of PNS convection. 
Second, PNS convection affects the evolution of the gain layer 
and shock due to a less compact PNS and thus larger neutri- 
nospheric radii with slightly lower temperatures and signifi¬ 
cantly higher luminosities near the neutrinospheres. This leads 
to a different structure between neutrinospheres and shock even 
without HB convection. 

3.2. Hot-bubble convection 

We now turn to the discussion of Ledoux-type convection in 
the gain layer behind the shock, the so-called “hot bubble” 
(HB) convection. Convective instability starts to develop in 
this layer in Models sl5_32 and s20_32 at about 30 ms after 
bounce (Fig. 8), and convective anisotropies become visible, 
e.g. in the velocity and entropy distributions, shortly afterwards 
(see Fig. 18). The radial structure of the star, the shock radius 
(Fig. 15), and thus the evolution of the models after bounce 
are, however, hardly affected. In this sense the convective ac¬ 
tivity in the hot-bubble region is “weak”. The “growth num¬ 
ber” Hgrow (Efi- 4, Fig. 3) for the growth of perturbations in 
the advection flow between shock and gain radius is corre¬ 
spondingly low. Obviously, neutrino heating is not powerful 
enough to drive strong buoyancy against the inward motion of 
the advected fluid. The small shock radius is associated with 
high negative postshock velocities, thus damping the growth of 
convection (for a general discussion, see Foglizzo et al. 2005). 
Only during the transient shock expansion, which sets in when 
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Fig. 17. First (top) panel: Mass in the gain layer. Second panel: 
Advection timescale as defined by Eq. (8). Third panel: Total heating 
rate in the gain layer. Fourth panel: Ratio of advection timescale to 
heating timescale. All lines are smoothed over time intervals of 5 ms. 
Note that the evaluation of Tadv fails after the onset of the expansion 
of the postshock layer in the exploding Model sll2_128_f. 

the entropy jump at the Si-SiO interface reaches the shock, 
does the advection timescale become sufficiently long to give 
rise to a phase of slightly stronger overturn motions. 

More powerful HB convection occurs in Model sll2_64. 
As we can see from Fig. 3, the perturbation growth factor in¬ 
creases to values of more than five in simulations with the si 12 
progenitor when the Si-SiO interface passes the shock 90 ms 
after bounce. The subsequent expansion of the shock is convec- 
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Fig. 18. Postshock convection in Model sl5_32. Figure a shows snap¬ 
shots of the entropy (in k^Poy) at four post-bounce times. Figure b 
shows the radial velocity at the same times (the color bar is in units of 
1000 km s"^ and its upper end corresponds to 2500 km s"^). The equa¬ 
torial plane of the polar grid is marked by the black diagonal lines. 


Fig. 19. Postshock convection in Model si 12_64. Figure a shows snap¬ 
shots of the entropy (in k^jhy) at four post-bounce times. Figure b 
shows the radial velocity (in 1000 km s"^) at the same times, with 
maximum values of up to 35000 km s"^ (bright yellow). The equato¬ 
rial plane of the polar grid is marked by the diagonal lines. 


lively supported and much stronger than in the ID model (see 
Fig. 15). Although there is initially no obvious morphologi¬ 
cal difference in the convective pattern before and right after 
the shock expansion (compare Fig. 19, upper left panels with 
upper right panels), the conditions in case of a larger shock 
radius become more favorable for a significant strengthening 
of the neutrino-driven convection, in contrast to the situation 
in the more massive progenitors, where the composition inter¬ 
face has a much smaller impact. The long advection timescale 
of 30-40 ms (Fig. 17, 2nd panel) leads to a ratio of advection 
to heating timescale close to unity (Fig. 17, fourth panel). The 


mass in the gain layer with specific energy above certain val¬ 
ues increases (Fig. 20), the pressure increases, and the shock 
is pushed to larger radii. The shock is highly deformed and the 
average shock position therefore fluctuates strongly due to the 
vigorous convective activity. 

We find that in particular large-scale modes of the flow pat¬ 
tern gain strength when the gain layer becomes radially more 
extended. Besides Ledoux-instability of the neutrino-heated 
layer, the growth of these modes can be supported by the 
vortical-acoustic cycle in accretion flows (cf. Foglizzo 2001, 
2002 and Foglizzo et al. 2005), which has recently been dis- 
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Fig. 20. Mass in the gain layer with the local specific binding energy 
as defined in Eq. (3) (but normalized per nucleon instead of per unit of 
mass) above certain values. The results are shown for Models si 12_32 
and sll2_64 for times later than 100 ms after core bounce. 


covered at work during the accretion phase of the stalled super¬ 
nova shock (Blondin et al. 2003; Scheck et al. 2004; Blondin & 
Mezzacappa 2006; Ohnishi et al. 2006; Scheck et al., in prepa¬ 
ration). The typical wavelength of the most unstable convective 
modes turns out to be linked to the radial width of the con¬ 
vective shell (for details, see Foglizzo et al. 2005). In Model 
sll2_64 the number of rising high-entropy plumes decreases 
from five at 50 ms to only two at 107 ms and even only one 
at 115 ms after bounce (the equatorial bubble in Fig. 19 at 
107 ms merges with the other bubble shortly afterwards). Such 
big structures persist for about 20 ms before they collapse and 
new rising bubbles form. The pattern is very nonstationary. 
Also at 140 ms and 170 ms large single bubbles appear, which 
dissolve again after roughly 20 ms. After the third generation 
of large bubbles has disappeared, the outer radius of the con¬ 
vective layer shrinks because of shock contraction, and smaller 
convective cells form again. 

The sequence of generations of big, floating bubbles pro¬ 
duces an oscillatory behavior of the shock radius (Fig. 15) 
as a consequence of a kind of feedback cycle between neu¬ 
trino heating, bubble expansion, shock expansion and over¬ 
shooting, reduced heating, shock contraction, bubble compres¬ 
sion and collapse, increased accretion and heating, new bubble 
growth, and again shock expansion. This is similar to the os¬ 
cillatory phases of shock expansion and recession, which we 
observed for one of the ID models discussed in Buras et al. 
(2006, Sect. 3.1.4). In the present 2D situation, however, there 
is not the same coupling between shock behavior and neu¬ 
trino luminosities found in that ID model, where the accretion 
flow and thus the neutrino luminosities and neutrino heating 
were quenched during phases of rapid shock expansion. In the 
multi-dimensional case, downflows of accreted matter around 


a rising plume can still feed the cooling layer near the neu- 
trinosphere and keep the neutrino luminosities at a high level 
(see Fig. 12). Here the expansion of the dominant rising plume 
lowers the heating rate of the bubble material, because its gas 
moves quickly away from the region of strongest neutrino en¬ 
ergy deposition. This means that the expanding plume cuts off 
its own energy supply. The expansion is, however, still too 
weak to push the shock outward against the ram pressure of 
the infalling gas of the progenitor star. After transiently over¬ 
shooting a possible equilibrium position, the shock therefore 
turns around, causing the compression of the high-entropy bub¬ 
ble and thus the onset of rapid cooling of the hot bubble gas 
by neutrino emission. The bubble collapses, and cool gas is 
channelled from the shock to the gain radius, increasing again 
the mass in the region of strongest heating. The energy input 
produces a new rising plume, which drives the shock outward 
again. The crucial question is whether the energy in the gain 
layer and the amount of matter with high specific energy in¬ 
creases or decreases during such a feedback cycle. In the first 
case an explosion may be produced as it happened in the ID 
model of Buras et al. (2006) (where it was favored by a ma¬ 
nipulation of the neutrino transport). In the present case, how¬ 
ever, the feedback cycle ends after three periods (these three 
periods are visible in Fig. 15 as well as in Fig. 20). The de¬ 
cisive quantity is the ratio of advection to heating timescale, 
the time-average of which should be larger than unity to obtain 
net energy gain in the heating layer. For Model si 12_64 this is 
not the case (Fig. 17, bottom panel). Although Model sll2_64 
gets much closer to success than any of the other simulations 
because of its strong convection and the large timescale ratio 
of almost unity, the simulation Anally still fails to produce an 
explosion. 

Turning now to resolution studies, we can compare Model 
sll2_64 with Model sll2_32. The latter has only half the lat¬ 
eral resolution but is otherwise identical with respect to input 
physics and radial grid. Clearly, the shock does not expand as 
far as in the high-resolution model (Fig. 15) and convection 
does not become as vigorous. Although the tendency of one 
convective cell to form with one rising plume and one down¬ 
flow can also be observed in the less resolved model, the con¬ 
vective activity appears weaker and less neutrino-heated matter 
is involved in the overturn motions (Fig. 20). One reason for 
this is the fact that downflows become too narrow to be well 
resolved in the low-resolution model. Their lateral width and 
stability are therefore numerically overestimated, thus chan¬ 
nelling more mass into the cooling layer below the gain radius. 
Numerical viscosity has a visible influence on the calculation. 

Although the resolution in Model si 12_64 is probably still 
not sufficient, the crucial quantity, the timescale ratio, increases 
only insignificantly when the lateral resolution is improved 
from 2.7 degrees to 1.35 degrees. We therefore have doubts 
that a further increase of the resolution can bring the timescale 
ratio above unity and thus establish favorable conditions for an 
explosion. Anyway, the resolution in Model sll2_64 appears 
to be the minimal resolution needed for reasonably converged 
simulations. An angular cell width of one degree would be 
preferable according to systematic resolution tests (L. Scheck, 
private communication). 



18 


R. Buras et al.: Two-dimensional hydrodynamic core-collapse supernova simulations with spectral neutrino transport 




25 

20 

15 

10 

5 




25 

20 

15 

10 

5 


a 


-200 


0 

r [km] 


200 








cn 




cn 


c 



r [km] 



Fig. 21. Postshock convection in Model sll2_128_f. Figures a-c show snapshots of the entropy for six post-bounce times. Figures d-f display 
the radial velocity at the same times with maximum values of up to 47000 km s"^ (bright yellow). In Figure d also the convective activity below 
the neutrinosphere (at radii r ^ 30 km) is visible. It is characterized by small cells, which are very similar to those observed in all other 2D 
models, too. The polar axis of the spherical coordinate grid is directed horizontally, the “north pole” is on the right side. 















R. Buras et aL: Two-dimensional hydrodynamic core-collapse supernova simulations with spectral neutrino transport 


19 


The size of the seed perturbations for starting convection 
turns out not to have an important influence on the HB con¬ 
vection. Comparing Models sl5_32 and sl5_64_p we And that 
a larger seed (as in the latter model) leads to HB convection 
which reaches the nonlinear regime about 10 ms earlier. The 
subsequent evolution, however, is very similar, and the faster 
onset of convective activity does not have any noticeable long¬ 
time effects. The reason may be the fact that the other rele¬ 
vant timescales, in particular Theat and Tadv, are typically longer 
than the perturbation-dependent differences of the growth of 
HB convection. 


3.3. A simulation with a grid of 180 degrees 

We now discuss Model si 12_128_f. The model is different from 
Model sll2_64 in some important technical aspects. The full 
star was simulated with a lateral grid of 180° instead of the 
~ 90° wedge around the equatorial plane, which we used for 
Model sll2_64. Along the polar axis of the spherical grid the 
gas flow is reflected, in contrast to our choice of periodic con¬ 
ditions at the lateral boundaries of the 90° wedge. Due to the 
assumed symmetry in case of the 90° grid, artificial constraints 
are imposed on the fluid flow. This has a selective influence 
on the type of flow pattern which can develop in the neutrino- 
heated hot-bubble layer. Global asymmetries with a dominant 
contribution of low modes — i.e., flow with 1=1 (dipolar) or 
I = 2 (quadrupolar) character in terms of orders of an expan¬ 
sion in Legendre polynomials for the cosine of the polar angle 
— can, for example, only be seen when the star as a whole 
is simulated. The possibility of long-wavelength modes in the 
postshock convective zone of the 11.2 M© simulation is in fact 
suggested by the linear perturbation analysis of Foglizzo et al. 
(2005), because ^grow of Eq. (4) reaches values around 6-7 in 
our ID run (Fig. 3). Low-mode convection and global asym¬ 
metries of the flow morphology are indeed found to grow and 
to cause significant quantitative differences in the evolution of 
Model sll2_128_f compared to Model sll2_64. These differ¬ 
ences are large enough to change the outcome of the simulation 
for the 11.2 M© progenitor even qualitatively. 

Noticeable differences between both models occur after 
about 60 ms post-bounce evolution (see Fig. 21). Although the 
overall morphology of the flow pattern in the convective post¬ 
shock layer looks still similar in both simulations, bubbles at 
the poles of Model sll2_128_f extend to larger radii. This is 
probably a consequence of the numerical setup of the simula¬ 
tion, which assumes axial symmetry and a polar grid with an 
axis that is impenetrable for the fluid flow. Flow which con¬ 
verges towards the axis is directed either inward or outward. 
Therefore its motion and behavior are constrained by the ex¬ 
istence of the polar grid axis. Moreover, polar features have 
smaller volumes compared to structures near the equator, which 
are treated as tori around the symmetry axis. This geometrical 
difference of polar and equatorial structures is known to lead 
to differences in the growth rate of perturbations. Laser exper¬ 
iments (Kane et al. 2000) as well as 3D simulations (Scheck 
2005) suggest that the smaller axial bubbles can grow faster 
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Fig. 22. Shock radii at the poles and at the equator versus post-bounce 
time for Model si 12_128_f. 

than equatorial tori and behave more like the mushrooms of 
Rayleigh-Taylor instabilities in the truely 3D case. 

In the further evolution of Model sll2_128_f the shock 
and convective layer develop a large deformation with domi¬ 
nant dipole and quadrupole modes, which become more and 
more prominent. Huge bubbles inflate alternatingly in both 
hemispheres, while downflows are present near the equatorial 
plane. These downflows are very nonstationary and flutter back 
and forth between the hemispheres. Eventually, 180 ms after 
bounce, the gain layer is completely dominated by the two po¬ 
lar bubbles and one or sometimes two downflows around the 
equator, separated by a small transient bubble in between. The 
shock in this model reaches a radius of about 600 km at the end 
of our simulation (Figs. 15 and 22) and is expanding then with 
a speed that is typical of exploding models (~ 10000 km s“^). 
Unfortunately we had to terminate the simulation at 225 ms af¬ 
ter bounce before it was possible to deduce the final parameters 
of the beginning explosion. The simulation had to be stopped 
because of a lack of computer time and the small timesteps en¬ 
forced by large and rapid fluctuations of physical variables in 
the region where the equatorial downflow penetrates deep into 
the neutrino spheric layer. 

The polar plumes can be weakened transiently because of 
the ram pressure of the infalling material ahead of the shock 
or because of a descreasing supply of neutrino-heated matter 
when the downflows feed the high-entropy lobe in the oppo¬ 
site hemisphere. But the polar bubbles never collapse and con¬ 
traction phases are reversed by new, powerful waves of high- 
entropy, high-velocity matter expanding away from the equa¬ 
torial plane. The snapshot at 225 ms in Fig. 21 shows such a 
phase for the hemisphere on the left side (the polar grid axis is 
oriented horizontally): Above 400 km the matter is still expand¬ 
ing, but below it has started retreating due to the ram pressure 
exerted by the matter falling through the shock. However, from 
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Fig. 23. Kinetic energy, as a function of post-bounce time, associated 
with the lateral velocities of the matter between gain radius and shock 
(solid) and between the Ve sphere and shock (dashed) for Models 
si 12_64 and si 12_128_f. The rapid increase in case of Model si 12_64 
at tpb ^ 180 ms is a numerical effect associated with lateral fluid flows 
which are enabled by our use of periodic boundary conditions at the 
angular grid boundaries. 
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Fig. 24. The upper panel displays the “explosion energy” of Model 
sll2_128_f, defined by the volume integral of the “local specific bind¬ 
ing energy” as given in Eq. (3), integrated over all zones where 
this energy is positive, i.e. Fexpi = f ^ 0)pdV. The lower panel 
gives the mass in the gain layer with the local specific binding energy 
(per nucleon) above certain values. The evolution of these quantities 
as functions of post-bounce time is shown. Note the difference from 
Fig. 20. 


below 300 km a new plume with very high velocities is already 
reviving the expansion of the gas between 300 km and 400 km. 
The whole situation is very dynamical. 


The dipolar expansion is therefore driven and powered by 
the flow of neutrino-heated gas that is continuously replenished 
near the gain radius by the equatorial downflow of accreted 
matter. The consequences for the dynamical and energetic evo¬ 
lution of Model si 12_128_f can be easily inferred from inspect¬ 
ing global quantities. In Fig. 17, bottom panel, we see that 
the timescale ratio T*^^/Theat reaches values above unity in this 
model. A crucial difference of Model sll2_128_f compared to 
its 90° counterparts is the fact that the advection timescale does 
not decline again after its maximum as in the other models but 
remains essentially constant for t ^ 100 ms post bounce. After 
a period of the order of ~ r^dv the timescale ratio starts ris¬ 
ing monotonically to climb to nearly 1.4 until 180 ms, associ¬ 
ated with an expansion of the average shock radius (Fig. 22). 
This means that the dipolar instability has a stronger influence 
than Ledoux convection in the other models and leads to an en¬ 
hancement of the efficiency of neutrino energy deposition. The 
total energy transfer by neutrinos becomes sufficiently large (in 
fact rises at r ^ 170 ms after bounce; Fig. 17) so that some part 
of the matter in the gain layer becomes nominally unbound and 
strong expansion sets in. Figure 24 shows that the mass of the 
gas with total local specific energy above some limits increases 
continuously, and the mass of the matter with positive energy 
follows this trend. Therefore we detect a steep growth of the 
“explosion energy” (Fig. 24, upper panel) after 180 ms. This 
moment coincides with the accelerated expansion and marks 
the onset of the explosion. In Fig. 12 we see that the neutrino 
luminosities do not decrease compared to the luminosities of 
the non-exploding Model sll2_64. This suggests that in spite 
of the launch of the explosion neutrino heating will go on at 
a significant level and will deliver more energy to the ejecta. 
Therefore Model sll2_128_f does not exhibit the disadvanta¬ 
geous situation of explosions in spherical symmetry where the 
start of rapid shock expansion quenches the accretion of the 
forming neutron star and thus leads to a significant reduction 
of the neutrino luminosity. 

The development of large-scale anisotropies with I = 
1,2 modes was also seen in other 2D simulations (Blondin 
et al. 2003; Scheck et al. 2004; Ohnishi et al. 2006; Burrows 
et al. 2006) and occurred in some preliminary 3D simulations 
(Mezzacappa & Blondin 2003, L. Scheck, private communi¬ 
cation), too. The morphology of our 180° model looks similar 
to the models published by Blondin et al. (2003), who stud¬ 
ied the hydrodynamics of nonradial instabilities of accretion 
shocks without taking into account the effects of neutrino trans¬ 
port and neutrino heating. However, we do not find that the 
growth of “turbulent energy” in the expanding layer behind 
the shock plays a significant role for the explosion of Model 
sll2_128_f. In contrast to the simulations by Blondin et al. 
(2003), Model si 12_128_f shows a saturation of the kinetic en¬ 
ergy of the lateral gas motions between neutron star (or gain 
radius) and shock at times later than ~ 140 ms after bounce on 
roughly the same level"^ as in Model sll2_64. We neither find 
specifically large values of £’kin,iat in case of the 180° simula- 


^ The large values of Fkinjat between the neutrinosphere and the 
gain radius in case of Model sll2_64, which can be inferred from 
Fig. 23, are a consequence of our choice of periodic boundary condi- 
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tion, nor do we see a continuous growth or a distinct increase 
of that quantity correlated with the onset of the explosion and 
the development of positive values for the explosion energy af¬ 
ter 180 ms (Fig. 23). We therefore conclude that the explosion 
of Model sll2_128_f is driven by neutrino energy deposition 
and not by the amplification of turbulent kinetic energy in the 
postshock layer. 

Dipolar asymmetries seem to be generic phenomenon in su¬ 
pernova simulations in which the neutrino heating is not strong 
enough to revive the stalled shock within a short time, i.e. on 
a timescale of the order of the advection timescale in the gain 
layer. However, the question has to be asked whether the mor¬ 
phology of the fiow is enforced in the 2D simulations by the 
choice of the azimuthal symmetry and the coordinate grid with 
its polar axis. The existence of a preferred axis direction pre¬ 
determines and might even support the growth of a dipolar 
asymmetry in this direction. In contrast, 3D simulations with¬ 
out a preferred axis might reveal a slower development of such 
asymmetries. It therefore needs to be demonstrated that free¬ 
dom of motion in all directions does not prevent the instabili¬ 
ties with dominant I = 1, m = 0 contributions. In reply to such 
objections one might say that in case of only a small amount 
of rotation the polar axis will obtain a physical meaning and 
will act as a centrifugal barrier for the fiow (if the axial compo¬ 
nent of the specific angular momentum of fiuid elements is con¬ 
served). The direction for the development of dipolar deforma¬ 
tion might therefore naturally be selected by the axis of (slow) 
rotation. More arguments in favor of the physical nature of low¬ 
mode, global asymmetries in the supernova core come from an¬ 
alytic investigations and linear analyses. On the one hand these 
suggest that stalled accretion shocks are subject to nonradial 
instabilities with highest growth rates for / = 1 deformation 
(Thompson 2000; Foglizzo 2002). On the other hand it was 
found that the / = 1 mode is the fastest growing mode in case of 
volume-filling thermal convection in a sphere (Chandrasekhar 
1961). The latter situation might apply to the supernova core, 
provided the shock radius is much larger than the inner bound¬ 
ary of the neutrino-heated, convective layer. But the true su¬ 
pernova conditions differ in important aspects, e.g., boundary 
conditions and accretion fiow in the postshock layer, which pre¬ 
vent a direct quantitative application of Chandrasekhar’s anal¬ 
ysis (Foglizzo et al. 2005). In any case, 2D simulations as well 
as analytic and linear analysis can only be taken as suggestive. 
The relevance of low-mode asymmetries during the phase of 
shock revival will ultimately have to be demonstrated by 3D 
models without the handicap of a coordinate singularity along 
the polar axis of a spherical or cylindrical grid. 

Therefore, we have to cautiously interpret the results of our 
present simulation. The choice of the lateral grid size and of the 
boundary conditions has turned out to decide whether the dipo¬ 
lar instability is suppressed as in case of the 90° wedge with 
periodic boundary conditions, or whether it can develop (and 
possibly is supported or even enforced) as in case of the 180° 
simulation with axial symmetry and reflecting boundary con¬ 
ditions along the grid axis. While Model 112_64 with the 90° 


tions at both lateral boundaries of the 90° wedge used in this model, 
see Seet. 3.1. 


wedge fails marginally. Model si 12_128_f yields an explosion 
for the considered 11.2 Mq progenitor star. Because of the rel¬ 
atively low mass in the gain layer, the explosion might remain 
rather weak, but our simulation had to be terminated too early 
to be finally conclusive in this point. The quantitative similar¬ 
ity in many aspects but qualitative difference in the outcome 
of the two simulations demonstrates how close Model 112_64 
was already to an explosion. Conversely, it can also mean that 
small effects which weaken the growth of the dipolar instability 
might delay the onset of the explosion (and thus change quan¬ 
tities such as the initial mass cut and the explosion energy), or 
might even lead to a failure^. 

The bottomline is that our 11.2 M© simulation lingers at the 
borderline between failure and explosion. But the effect which 
has triggered the success in case of Model sll2_128_f, i.e. the 
large-scale non-radial modes of the fiuid fiow between shock 
and gain radius, can be treated only approximately in two- 
dimensional simulations. We can not exclude that the success 
is a result of an overestimation of this phenomenon and of its 
consequences. Three-dimensional (axis-free) simulations with 
reliable neutrino transport are needed to convincingly demon¬ 
strate that the neutrino-driven mechanism, supported by low¬ 
mode convection and accretion shock instability, is viable to 
explain supernova explosions of massive stars at least for some 
range of progenitor masses. 

3.4. A simulation with rotation 

We investigate the effects of rotation on the supernova evolu¬ 
tion with Model sl5_64_r, for which we assume that the iron 
core rotates with a period of about 12 s prior to collapse. This 
is rather slow compared to the conditions considered in many of 
the core-collapse simulations performed by other groups (e.g., 
Ott et al. 2004; Kotake et al. 2004). The reasons for our choice 
of the initial rotation law were explained in detail at the be¬ 
ginning of Sect. 3 and in Muller et al. (2004). The angular fre¬ 
quency Q of the Fe-core at the beginning is 0.5 rad s“^ (Fig. 25) 
and increases due to angular momentum conservation during 
collapse to maximum values around 600-700 rad s“^ in the ho¬ 
mologous core shortly after bounce, see Figs. 25, 26a. Rotation 
generates a centrifugal force Fcent which at this time is at most 
6 X 10“^ of the gravitational force (Figs. 25, 26a). The 
ratio of rotational to gravitational energy, jS^ot = ^^rot/l^^gravL 
grows from an initial value of less than 10“^ to roughly 0.4% 
at bounce (Fig. 27). Because of the subsequent contraction of 
the PNS and the inflow of material with higher specific angu¬ 
lar momentum (y^), the ratio Fcent/Fgrav rises up to 25% near 
the equator (Fig. 25) and can reach values around 3 in the po¬ 
lar regions at a post-bounce time of 271.1 ms (Fig. 26b). The 
angular frequency increases up to maximum equatorial values 

^ In this context it is interesting to note that for an 11 M© progeni¬ 
tor from Woosley et al. (2002), which has a less steep density profile 
exterior to the iron core, Burrows et al. (2006) did not find an explo¬ 
sion before about 0.5 seconds after bounce. Besides a different pro¬ 
genitor, there were, however, also other differences. They, e.g., per¬ 
formed purely Newtonian simulations and used a flux-limited diffu¬ 
sion treatment for the neutrino transport, in which energy bin coupling 
by Doppler shifts and neutrino-electron scattering were ignored. 
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Fig. 25. Radial profiles of the specific angular momentum angular velocity and ratio of centrifugal to gravitational force, Fcent/Fgrav, in 
the equatorial plane versus (a) radius and (b) mass for Model sl5_64_r at different times. Note that the “enclosed mass” in multi-dimensional 
simulations is defined here as the mass within a spherical volume of specified radius. Quantities plotted versus enclosed mass therefore do not 
represent Lagrangian information. Time is normalized to bounce. Note that in the upper right panel the lines for times -175 ms, 0, and 10 ms lie 
on top of each other because of specific angular momentum conservation and negligibly small rotational deformation of the collapsing stellar 
core until later after bounce. 


around 2000 rad s ^, and even much larger values close to the 
polar axis (Fig. 26b). 

As a consequence, the forming neutron star develops an in¬ 
creasing degree of rotational flattening. The oblateness of the 
neutrinospheres is most pronounced during the later stages of 
the simulated post-bounce evolution (Fig. 26, panels c,d,f and 
Fig. 29), whereas the deeper interior of the PNS shows a much 
smaller deformation, and the convective layer at r ^ 25 km 
remains nearly spherical (Fig. 26, panels e,f). At the end of 
our simulation the polar value of the radius of the Vg-sphere 
is about 30 km, which is roughly the same as in non-rotating 
models, whereas the Vg-sphere is at 40 km near the equatorial 
plane (Figs. 11, 26). Also the rotation-generated differences of 
the gain radius in different directions grow with time, although 


the contrast visible in Fig. 16 shows the combined effects of 
rotation and local variations due to the convective anisotropies 
in the postshock region (see Fig. 26). The angular differences 
of the gain radius become roughly a factor of two until the end 
of our simulation for Model sl5_64_r and are thus somewhat 
larger than in the other models with strong convective over¬ 
turn in the neutrino-heating layer (in particular Models si 12_64 
and sll2_128_0. The radial structure of the rotating PNS dif¬ 
fers most strongly from that of the corresponding non-rotating 
Model sl5_32 near the equatorial plane, where the centrifugal 
force, Fcent, points in the opposite direction of the gravitational 
force, Fgrav Differences in the radial structure between the ro¬ 
tating and non-rotating PNSs are less significant at the poles. 
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Fig. 26. Snapshots of Model sl5_64_r. The rotation axis is oriented vertically. Panels a, b show the distributions of the specific angular momen¬ 
tum jz, of the angular frequency Q, and of the ratio of centrifugal to gravitational force for the moment of bounce and for a post-bounce time of 
271.1 ms, respectively. Note the different scales which represent enclosed mass M (which corresponds to values of mass enclosed by spheres 
of chosen radii) and radius r. The black lines mark the shock. Panels c, d depict quantities which are of interest for discussing the physics in 
the gain layer at the time of maximum shock expansion and at the end of the simulation, respectively. The thick black line marks the shock, 
the thin black lines indicate the gain radius and Ve sphere. Panels e, f provide analogous information inside the PNS at the onset of convection 
and at the end of the simulation, respectively. Instead of the gas entropy and Fe. the total entropy including neutrino contributions and the total 
lepton number are plotted. The three black lines in panel f indicate from outside inward the neutrinospheres of Ve, Ve. and v^. 
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Because of the rotational flattening the PNS as a whole is 
less compact than in the non-rotating case, i.e., the volume en¬ 
closed by the neutrinospheres is larger. The implications of a 
bigger average neutrinospheric radius for the global properties 
of the neutrino emission were discussed in detail in Sect. 3.1. 
Due to the lower temperatures at the more extended neutri- 
nosphere, the mean energies of the radiated neutrinos are re¬ 
duced compared to the corresponding ID model, associated 
with a decrease of the neutrino luminosities. In non-rotating 
models, this reduction holds only for the first ~ 150 ms after 
bounce and is compensated later. In fact it is over-compensated 
by the convective energy transport to the neutrino sphere, which 
leads to higher luminosities (in particular for muon and tau neu¬ 
trinos) in 2D simulations at r ^ 150 ms after bounce (in spite of 
continuously lower mean neutrino energies; Figs. 12 and 13). 
In the rotating Model sl5_64_r this “cooling effect” of the more 
extended neutrinosphere is significantly stronger than in the 
corresponding non-rotating Model sl5_32, causing even lower 
energies of the radiated neutrinos (Fig. 13). In fact, the neu¬ 
trinospheric temperatures are reduced so much by rotational 
expansion that this effect cannot be overridden by the convec¬ 
tive transport of energy to the neutrinosphere, in particular also 
because PNS convection tends to be slightly weaker in Model 
sl5_64_r than in Model sl5_32 (compare Fig. 26 with Fig. 5 
and see the discussion following later). The neutrino luminosi¬ 
ties in the rotating model stay therefore clearly below those of 
its non-rotating counterpart during ah of the simulated evolu¬ 
tion (Fig. 12), and the integral of the radiated energy is smaller 
than in the ID case (Fig. 14, lower panel), signalling that en¬ 
ergy is stored in rotation instead of being converted to neutrino 
emission. In contrast, the rotating model does not show a par¬ 
ticular behavior with respect to the lepton number loss (Fig. 14, 
upper panel), because the electron neutrino and antineutrino lu¬ 
minosities are both reduced by rotational effects in very similar 
ways (Fig. 12). The emission anisotropies caused by the ro¬ 
tational deformation of the neutrinosphere and by convective 
effects will be addressed in Sect. 3.5. 

In contrast to the oblate deformation of the PNS, the shock 
reveals a prolate shape during most of the computed evolu¬ 
tion, with a radius typically 50 km larger at the pole than at 
the equator (Figs. 15, 26, panels c,d). The polar bulge of the 
shock is maintained by a big convective bubble which exists 
rather stably between polar angles of about 10° and about 45°. 
After an initial, transient phase in which the convective activity 
grows and smaller convective cells merge to larger, volume- 
flhing structures, this large pole-near bubble comes to exist 
besides one or two additional plumes closer to the equator. 
Although these bubbles are again strongly time-dependent as 
in the other 2D models, and phases of bubble contraction are 
followed by bubble reinflation, the morphology is very stable 
and the bubbles grow again essentially at the same places. The 
pattern seems to be determined and supported by the presence 
of an angular momentum gradient in the neutrino-heated layer 
(Fig. 26, panel b) and by the action of centrifugal and coriolis 
forces on the fluid motion. 

Towards the end of our simulation (from about 180 ms un¬ 
til 280 ms post bounce) quasi-periodic large-amplitude pole- 
to-equator oscillations with a cycle time of 15-25 ms set in 
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Fig. 27. Ratio of rotational energy to the gravitational energy, y^rot. as 
a function of time for Model sl5_64_r. The quantity is shown for the 
whole computational volume (solid line) and for spherical volumes 
with an enclosed mass of the Fe-Si core (1.43 Mq, dashed line) and 
of the Fe core (1.28 Mq, dotted line), respectively. Note that in the 
latter two cases the time evolution of y^rot docs not provide Lagrangian 
information, because matter is redistributed by rotational deformation 
and by convection. 

(Fig. 15) where phases with a larger pole-near bubble alter¬ 
nate with time intervals in which the convective plume near 
the equator is stronger. The transition between both extrema 
is characterized by a merging of the two convective cells into 
one big plume at intermediate latitudes. The shape of the shock 
changes back and forth between a pronounced prolate deforma¬ 
tion and a more oblate shape. During ah these time-dependent 
variations the maximum and minimum shock radii stay around 
200 km and 150 km, respectively (Fig. 15). The average shock 
radius at r ^ 150 ms p.b. is 50-100% larger than in the non¬ 
rotating 2D models and does not decay until the end of our sim¬ 
ulation at nearly 300 ms after bounce. Due to centrifugal effects 
the presence of angular momentum in the infahing matter has 
a stabilizing influence on the postshock flow and on the shock. 
Thus rotation ensures a more extended gain layer and supports 
strong convection, in contrast to the non-rotating 15 M© Models 
sl5_32 and sl5_64_p, where the retraction of the shock after 
its maximum expansion strongly damps the convective activ¬ 
ity in the gain layer (Fig. 15 and Buras et al. 2006). The con¬ 
vective pattern in the postshock region therefore depends on 
the amount of angular momentum carried by the accreted mat¬ 
ter. Since the gas falling onto the shock comes from larger and 
larger initial radii at later times, the angular momentum of the 
accretion flow increases continuously (Fig. 25). This may ex¬ 
plain some of the evolution which we observe for the morphol¬ 
ogy and dynamical behavior of the convective post-shock layer. 

The conditions for convective instability are affected by ro¬ 
tation in a way which is expressed by the so-called Solberg- 
Hpiland criteria (see e.g. Tassoul 1978; Keil 1997). For the 
considered situation in the supernova core, the first criterion is 
a generalization of our Quasi-Ledoux criterion, combined with 
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the Solberg-term which accounts for the stabilizing effect of a 
positive angular momentum gradient. With the parameter yS^ifF 
for neutrino diffusion in the Quasi-Ledoux criterion (Eq. 7), the 
corresponding mode frequencies are 
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where instability is given when ojuoIQL > 0. The sign of this 
expression is identical to the sign of the absolute value term 
under the square root, ^cent is the acceleration due to cen¬ 
trifugal forces, Vtot is the sum of gas and neutrino entropies, 
and X = rsin?7 is the distance from the rotation axis. Here, 
the Quasi-Ledoux criterion was generalized further to a fully 
multi-dimensional form by introducing the scalar product of 
the vector quantities. In the gain layer, the adequate Solberg- 
Hpiland criterion is obtained from Eq. (9) by ignoring the neu¬ 
trino entropy and lepton number, i.e. by replacing Vtot by v and 
^lep by Yq, and setting = 0 when neutrino diffusion effects 
are irrelevant. The second Solberg-Hpiland criterion for insta¬ 
bility is given by 
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where z = r cos d is the distance from the equatorial plane and 
the suffix z denotes the z-component of vector quantities. We 
will not further consider this criterion in the following because 
convective instability requires at least one of the two criteria 
to be fulfilled. An evaluation of our model, however, reveals 
that the second criterion indicates instability only when the first 
criterion is also fulfilled. 

As already discussed in earlier publications (Keil 1997; 
Janka & Keil 1998; Janka et al. 2001) the differential rotation, 
which is accounted for by the first term in Eq. (9), tends to 
damp convective activity near the polar axis in the PNS, see 
Eig. 26e,f. The same effect is found for the hot-bubble layer, see 
Eig. 26c,d. At larger distances from the polar axis the distribu¬ 
tion of the specific angular momentum is flatter and convection 
is only weakly affected by rotation, i.e. convection develops in 
regions inside the PNS and between the gain radius and the 
shock basically similar to the situation in Model sl5_32. PNS 
convection therefore has effects analogous to those discussed 
in Section 3.1. The damping of convection near the poles leads 
to a slightly slower effective transport of lepton number (and 
energy), see Pig. 28 for the situation after 200 ms post-bounce 
evolution, especially at later times when the convection near the 
axis is more strongly suppressed due to a steeper gradient of 7 ^, 
see Pig. 26e,f. Another new feature in the rotating model is the 
transport and redistribution of angular momentum by convec¬ 
tion. Pig. 25b shows how convection in the PNS produces a flat 
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Fig. 28. Tiep profile at 200 ms post-bounce time for a ID model, a 
2D model, and our 2D model with rotation. The plot shows angle- 
averaged quantities for the 2D cases. 


radial profile of jz near the equator. This happens despite the 
stabilizing effect of the initially positive derivative of jz because 
the negative entropy gradient dominates the Solberg-Hpiland 
criterion and drives convective instability. 

How does rotation affect the possibility for getting explo¬ 
sions by the delayed neutrino-heating mechanism? The anal¬ 
ysis by Yamasaki & Yamada (2005) suggests that rotation can 
appreciably improve the conditions for shock revival in case the 
initial rotation frequency is at least / ~ 0.1 Hz at 1000 km. This 
corresponds to a rotation period of 10 s and is therefore only 
slightly faster than the rotation considered in Model sl5_64_r 
(where the initial spin period is 12 s in the iron core). Yamasaki 
& Yamada (2005) found a sizable reduction by 25% of the 
“critical neutrino luminosity” for starting a neutrino-driven ex¬ 
plosion at a mass accretion rate of about 1 M© s“^ or lower. The 
effects of rotation are, however, diverse and modify the struc¬ 
ture of the collapsing star, the convection in the core, the gas 
motion behind the shock, and the radius and neutrino emission 
of the forming neutron star. A separate discussion of selected 
effects as done by Yamasaki & Yamada (2005) can therefore be 
misleading. In our simulations all effects of rotation are fully 
coupled and we can assess the question how these effects in 
combination determine the conditions for the delayed explo¬ 
sion mechanism. To this end we again compare the results of 
our rotating model with the non-rotating counterparts. 

Rotation turns out to be helpful, but to a much lesser ex¬ 
tent than estimated by Yamasaki & Yamada (2005). Certainly, 
the shock radius is significantly larger than in the non-rotating 
models of the 15Mo star (Pig. 15), the effective advection 
timescale through the gain layer correspondingly becomes 
longer by up to a factor of about 4, and the mass in the gain 
layer, Mgi, increases to a value that is two or three times larger 
than in the non-rotating models. These differences are on the 
one hand an indirect consequence of the structural changes of 
the PNS (cf. Section 3.1), on the other hand they result directly 





26 


R. Buras et aL: Two-dimensional hydrodynamic core-collapse supernova simulations with spectral neutrino transport 



Fig. 29. a Mass Mns and total angular momentum /ns of the “neutron 
star” (defined by the matter with density above 10^^ g cm"^) in Model 
sl5_64_r. b Moment of inertia /ns and average angular frequency 
f^NS = ^NsZ-^NS of the neutron star in Model sl5_64_r. c Average ra¬ 
dius Rns of the neutron star, defined as the radius of a sphere with 
the same volume as the deformed rotating neutron star, “effective” ra¬ 
dius defined as RefF,NS = V^ns/^/ns, and neutron star radii at the pole 
(Rns,poi) and equator (RNS,eq). for Model sl5_64_r. Also shown is the 

eccentricity ^ns = -^1 “ ^Ns.poi/'^Ns.eq- 


from the influence of centrifugal forces on the fluid flow in the 


gain layer. In contrast, however, the total heating rate, and 
the ratio of advection to heating timescale increase only slightly 
(Fig. 17). The latter remains significantly below unity, because 
the heating timescale Theat has increased considerably due to the 
lower energies and luminosities of the neutrinos radiated from 
the rotating PNS (Figs. 12, 13). The rotating model therefore 
shows no tendency to develop an explosion until we stopped 
our simulation at about 280 ms after bounce. 

If angular momentum conservation holds during the subse¬ 
quent evolution, one can use our results at the end of the simu¬ 
lated evolution to roughly estimate the final angular frequency 
of the cold neutron star after its neutrino cooling and contrac¬ 
tion phase, QcNS, assuming it to be a rigid rotator. The “NS” 
is defined here as the mass at densities above 10^^ gcm“^. In 
Fig. 29 we show the time evolution of several quantities for 
this such defined “NS”, namely of the radius /?ns, defined as 
the radius of a sphere with the same volume as the deformed 
proto-neutron star, and of the average angular frequency, de¬ 
fined as Qns = ^ns/7ns, where /ns is the moment of inertia 
and /ns is the total angular momentum of the “NS”. Taking 
the conditions dX t = 270 ms after bounce and assuming an¬ 
gular momentum convervation — no processes happen which 
transport angular momentum out of the NS or add mass and 
angular momentum to it —, the angular frequency of the NS 
after self-similar contraction (i.e., the shape of the star does not 
change and thus its moment of inertia scales with the square of 
the average NS radius) to a final radius /?cns = 10-12 km is 

7^ _ 7^ //?NS,sim\^ 

IZcNS - i^NS,sim ^- 

\ /<cNS / 

( 27 kin 

-———^^1 ^ 5 ... 7 rad ms“^, (11) 

10-12 km/ ^ 

where nNS,sim and /^NS,sim are the average angular frequency 
and average radius of the “NS” at the end of the simulation, 
respectively (Fig. 29). This corresponds to a spin period of 
Fcns = 27r/QcNS - 1 ms, which is at least 10-100 times faster 
than the typical birth period of pulsars estimated from observa¬ 
tions (see Heger et al. 2005 and references therein). Therefore 
our Model sl5_64_r is on the extreme side of conditions which 
are present in collapsing stellar cores in typical cases (see also 
the discussion in Ott et al. 2006). 

3.5. Anisotropy of the neutrino emission 

In this section we shall discuss the variations of the neutrino 
emission with polar angle, which are a consequence of PNS 
convection, convective overturn in the layer between gain ra¬ 
dius and shock, and of rotation. For this purpose we will mostly 
concentrate on the non-rotating Model sll2_128_f, which de¬ 
velops extremely strong convection and anisotropy, and on 
our case with rotation. Model sl5_64_r. The anisotropy of the 
neutrino emission from rotating collapsing stellar cores was 
first evaluated by Janka & Monchmeyer (1989a,b) analytically 
and by means of post-processing Monte Carlo transport cal¬ 
culations using 2D core-collapse models. They found that the 
neutrino emission is higher in the polar direction than in the 
equatorial plane, verifying von ZeipeTs law of gravity dark- 
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Fig. 30. a Lateral maxima of the Ve flux relative to the average Ve flux 
versus time for Models sl5_64_r and sll2_128_f, evaluated at 400 km 
radius for an observer at rest, b Time integral of the “isotropic equiv¬ 
alent luminosity” Ljso - calculated by assuming that the flux in one 
polar direction is representative for all other directions - versus cos 
for different models and neutrinos. The integral was performed from 
10 ms to 200 ms after bounce; cos?7 = +1 correspond to the poles, 
cos = 0 to the equator. The thin lines are for an evaluation at 400 km 
radius for an observer at rest, the thick lines are given only for Ve and 
show the result at the Ve-sphere for an observer at rest. 


ening in case of neutrinos from supernova cores. This result 
was confirmed by Kotake et al. (2003), who also performed 
a post-processing analysis of 2D core-collapse calculations in 
which a trapping scheme for the neutrino treatment was em¬ 
ployed. It must be pointed out, however, that none of these ap¬ 
proaches was self-consistent: The emission anisotropy and the 
neutrino treatment in the hydrodynamic simulations were cal¬ 
culated with different approximative methods, the feedback of 
neutrino transport on the hydrodynamics was not taken into ac¬ 
count (a trapping scheme only computes local source terms), 
and the evaluation was based on simplifying concepts, e.g., 
the assumption that a well-defined average neutrinosphere ex¬ 



a cos'll 



Fig. 31. a Ratio of electron neutrino flux to average flux (for an ob¬ 
server at rest at 400 km) versus cosine of the polar angle for Model 
sll2_128_f for different post-bounce times. The times are picked such 
that large maxima of the flux ratio occur (see Fig. 30a). b Same as 
panel a, but at the Ve sphere. The vertical solid lines mark the lateral 
boundaries of the region where the results were evaluated for Fig. 30a. 
Higher and lower latitudes were excluded because of possible artifacts 
from the reflecting boundaries. 

ists^. The calculated anisotropies must therefore be interpreted 
with great caution and conclusions drawn on their basis may 
be very misleading, in particular concerning their implications 
for the explosion mechanism. The latter depends on many com¬ 
peting effects and a discussion requires a coupled and consis¬ 
tent treatment of neutrinos and hydrodynamics. Walder et al. 
(2005) recently performed such simulations of rotational core 
collapse, in which 2D flux-limited neutrino diffusion was fol¬ 
lowed self-consistently in 2D hydrodynamic models. These au- 

^ The Monte Carlo calculations by Janka & Monchmeyer (1989b) 
did not make use of such simplifications, assumed, however, that the 
transport in different lateral directions of the 2D environment can be 
approximated by a “ray-by-ray” approach, in which the neutrino flux 
at a given latitude is assumed to be radial and can be calculated by 
the transport in a spherical background with a radial structure as it is 
present in the 2D model at this latitude. 
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b cosi^ 


Fig. 32. Same as Fig. 31, but for Model sl5_64_r. The pole of the ro¬ 
tating model corresponds to cos = 1. 

thors pointed out the fact that the neutrino emission is always 
more isotropic than the rotationally deformed or convectively 
distorted mass distribution because of the non-locality of the 
transport, i.e., the neutrino distribution at every point is a su¬ 
perposition of the irradiation from different contributing direc¬ 
tions. The paper by Walder et al. (2005), however, does not 
provide quantitative information in a form which would allow 
us to make comparisons with their results. In the following we 
will therefore concentrate on the presentation and discussion of 
our models. 

Figure 30a shows the lateral maximum of the ratio of the lo¬ 
cal neutrino flux to the average flux for Vq as function of time. 
In case of Model sll2_128_f spikes are visible with a typical 
width of about 5 ms and an amplitude that grows with time 
and reaches maximum values up to 2.4 towards the end of our 
simulation. These spikes are caused by “outbursts” of neutrino 
radiation, which are associated with very narrow downflows 
of accreted matter, which reach deep into the cooling layer 
and create there local “hot spots” (Fig. 21). These downflows 
and their location, however, are extremely time-dependent and 
nonstationary. This can be seen in Fig. 31a where the lateral 


flux variations are shown for three different times after core 
bounce. The large-amplitude flux variations are mostly a phe¬ 
nomenon which affects the neutrino emission from the cooling 
layer between gain radius and neutrino sphere. The amplitudes 
of / (F) are much smaller at the neutrino sphere (compare 
Figs. 31a andb). 

Because the downflows are transient in space and time, the 
lateral variations disappear in the time-integrated flux at dif¬ 
ferent latitudes, and only minor fluctations remain (Fig. 30b). 
Local maxima or minima at the lateral edges are a consequence 
of the reflecting boundary conditions, which cannot be pene¬ 
trated but redirect the flow in the inward or outward direction. 

The flux ratio maxima for Model sl5_64_r also exhibit 
short-time variations with typical fluctuation timescales of 
about 5 ms, but with lower amplitudes (up to values of 1.5) than 
those in Model sll2_128_f. This is mainly due to the fact that 
the convective layer in Model sl5_64_r has a smaller radial ex¬ 
tension, the downflows are less narrow and hit the cooling layer 
with less violence, producing less extreme local emission than 
in Model si 12_128_f. An exception to this is the polar region of 
the rotating model, where the accretion flow is able to create a 
high-luminosity spot even at the neutrino sphere, causing larger 
latitudinal flux variations there than farther out (Figs. 32a,b). 

While the time-integrated flux at large distances reveals a 
local maximum near the equator in Model sl5_64_r (Fig. 30b), 
which is associated with a persistent equatorial downflow, the 
time-integrated flux at the neutrino sphere is nearly featureless 
with only a very shallow global pole-to-equator gradient. The 
latter is a consequence of the rotational flattening of the PNS, 
which, however, is too low to have significant effects on the 
instantaneous (Figs. 32b) or time-integrated (Fig. 30b) neutri- 
nospheric emission. 

We point out that a discussion of the lateral variation of the 
neutrino emission of our models is handicapped by the fact that 
our approximation of 2D neutrino transport disregards the lat¬ 
eral component of the neutrino flux vector and therefore tends 
to overestimate the angular asymmetry of neutrinos streaming 
out from radiating regions (for a more detailed discussion, see 
Buras et al. 2006). Truely multi-dimensional transport should 
therefore not only reveal smaller lateral variations of the time- 
integrated flux as our models do, but will also show less ex¬ 
treme angular variations of the instantaneous emission on short 
spatial wavelengths. 

4. Summary and conclusions 

We have presented results of a series of core-collapse and 
post-bounce simulations for different progenitor stars between 
11.2 Mo and 25 Mq, comparing 2D (axially symmetric) with 
ID (spherically symmetric) calculations. Doing so, our main 
goals were (i) investigating the differences between convection 
in progenitors with different masses, (ii) exploring the effects 
of convection below the neutrino sphere (“PNS convection”) on 
the proto-neutron star structure, its neutrino emission, and the 
neutrino heating-layer behind the shock, (iii) investigating the 
role of hydrodynamic instabilities that affect the stalled accre¬ 
tion shock, i.e. convective overturn in the neutrino-heated “hot 
bubble” layer (“HB convection”) and global low-mode nonra- 
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dial instability of the accretion shock (termed SASI by Blondin 
et al. 2003 and possibly caused by the action of an advective- 
acoustic cycle according to Foglizzo 2001, 2002), (iv) studying 
the effects of rotation, and (v) testing the influence of numeri¬ 
cal aspects like the grid resolution, size of the angular wedge, 
and magnitude of seed perturbations for convection. Since our 
2D neutrino-hydrodynamics code is a direct descendant of our 
ID Prometheus/Vertex code, it is particularly well suited for 
performing such comparisons of ID and 2D supernova models. 

Convection inside the proto-neutron star starts 30-40 ms af¬ 
ter bounce in all of our 2D models and encompasses a layer 
growing in mass until the end of our simulations (which were 
typically terminated about 250 ms after bounce). It leads to a 
more extended neutron star than in the ID simulations with 
lower temperatures at the neutrinosphere. For this reason the 
mean energies of the neutrinos emitted from the neutrinosphere 
are reduced (up to 10% after 200 ms of PNS convection). 
Despite the larger radiating surface, the lower neutrino spheric 
temperatures also cause a slight reduction of neutrino luminosi¬ 
ties during the first 150 ms after bounce. This holds in partic¬ 
ular for Ve, because the convective transport of lepton number 
maintains a higher electron degeneracy in the neutrino spheric 
region and accelerates the lepton number loss compared to ID 
simulations. Only at r ^ 150 ms after bounce, convectively en¬ 
hanced energy transport in the nascent neutron star also leads 
to increased energy loss, and the luminosities of heavy-lepton 
neutrinos become significantly (15%-20%) higher than in the 
spherical models. 

PNS convection of the kind found in our simulations leads 
to a slightly reduced total energy deposition in the gain layer 
mainly because of the lower average energies of the radiated Ve 
and Vq. Since the effects of convection below the neutrinosphere 
are hard to disentangle from those of hydrodynamic instabili¬ 
ties in the neutrino-heating layer behind the stalled shock, we 
developed a simple “mixing algorithm”. It allowed us to re¬ 
produce ah major effects of PNS convection in ID simulations 
and thus to separate them from the consequences of multi¬ 
dimensional fluid flow in the postshock layer and to arrive at 
the above conclusion. 

Convective overturn in the neutrino-heating layer remained 
rather weak in case of the 15 M© and 20 M© progenitors. The 
main reason for that is the rapid contraction of the accretion 
shock after its maximum expansion. This causes the gain layer 
to be very narrow and the infah velocities of the gas ahead and 
behind the shock to be very high. As a consequence, the ad- 
vection timescale of the gas through the gain layer is very short 
compared to the typical neutrino-heating timescale. Buoyancy 
forces therefore hardly achieve bubble rise in the flow of gas 
accreted from the shock to the gain radius. As suggested by 
Janka & Keil (1998) and Janka et al. (2001) and verified by 
Thompson et al. (2005), the ratio of the advection timescale 
to the neutrino-heating timescale, Tadv/^heat, turned out to be 
a useful diagnostic parameter to measure the proximity of a 
model to a neutrino-driven explosion. A necessary condition 
for an explosion is that the timescale ratio rises above unity 
for a time interval of at least the neutrino-heating timescale. In 
case of the 15 M© and 20 M© models, HB convection increases 
the heating rate and the timescale ratio to values only slightly 


larger than in the ID simulations, but still roughly a factor of 
two below the critical limit. We therefore found explosions of 
these stars neither in spherical symmetry (in agreement with 
Liebendorfer et al. 2002, Thompson et al. 2003, Sumiyoshi 
et al. 2005) nor in 2D. 

Also rotation did not change this negative outcome. We 
studied one 15 M© model with pre-collapse rigid iron core ro¬ 
tation of ~ 12 s period, which leads to a neutron star with a spin 
period of about 1 ms, if the angular momentum of the core af¬ 
ter collapse is conserved. Rotation of this size is probably on 
the extreme side of what can be expected for the cores of “nor¬ 
mal” supernovae, which are supposed to give birth to neutron 
stars with an initial period of 10 ms or more (see the discus¬ 
sions in Heger et al. 2005 and Ott et al. 2006). Our simula¬ 
tions reveal a number of important differences of the rotating 
model compared to its non-rotating counterparts. The proto¬ 
neutron star develops an eccentricity of more than 0.6 until 
we stopped the simulation at nearly 300 ms after bounce. At 
this time the luminosities of the radiated neutrinos are signifi¬ 
cantly smaller (10-20%) and their mean energies up to 2 MeV 
lower than in the non-rotating 2D model, because the equa- 
toriahy more extended neutrinosphere is significantly cooler 
and energy is stored in rotation instead of being released by 
neutrinos. Despite the clear oblateness of the proto-neutron 
star, its rotation-induced emission anisotropy is very small. 
Nevertheless, rotation has a favorable influence on the condi¬ 
tions and parameters which determine neutrino-driven explo¬ 
sions. Centrifugal forces stabilize the accretion shock at larger 
radii, increase the advection timescale of the postshock gas sig¬ 
nificantly, and thus allow for a layer of well developed, strong 
convective overturn activity behind the shock. Because more 
mass stays in the gain layer for a longer time, the total energy 
deposition rate behind the shock is higher at later post-bounce 
times 130 ms p.b.) than in the non-rotating models. In spite 
of these healthy effects, however, the timescale ratio Tadv/^heat 
remains still well below unity (Tadv/^heat ^0.6). 

Even without rotation postshock convection becomes vio¬ 
lent in case of the 11.2M© star. The shock in this model is 
able to stay longer at large radii than in the more massive stars. 
This is due to the fact that the rather low-mass progenitor has 
a steeper density decline at the transition to the Si-rO layer, 
which leads to a rapid decrease of the mass accretion rate of 
the shock at about 90 ms after bounce. This allows the shock to 
reexpand in adjusting to the situation of reduced ram pressure. 
The increased advection timescale gives convection the possi¬ 
bility to gain strength and thus to support the shock at a much 
larger radius than in the corresponding ID model. Also the to¬ 
tal neutrino heating rate behind the shock and the efflciency 
of net neutrino-energy transfer to the gas in the gain layer is 
higher by up to a factor of two. The timescale ratio Tadv/'^heat 
approaches unity and remains close to — but slightly below 
— this threshold until the end of our simulations. The 11.2 M© 
model computed with a 90° lateral wedge therefore lingers at 
the border to success. 

Such a situation is extremely sensitive to relatively little 
changes. We saw this when we repeated the simulation with 
a full 180° grid instead of using the wedge around the equa¬ 
tor. While PNS convection turned out not to depend on the 
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wedge size, convective activity in the neutrino-heating layer 
can change significantly when the available degrees of free¬ 
dom are not constrained by periodic boundary conditions of a 
90° equatorial wedge and therefore low-mode deformation of 
dipolar (/ = 1) and quadrupolar (/ = 2) character is allowed for. 
Convection becomes sufficiently strong so that the accretion 
shock continues to expand. This ensures that the effective ad- 
vection timescale does not decrease after it has reached its max¬ 
imum. Ait 140 ms after bounce, the timescale ratio Tadv/Theat 
then becomes larger than unity, thus further improving the con¬ 
ditions for efficient energy deposition by neutrinos in the post¬ 
shock layer. After about 180 ms of post-bounce evolution the 
total energy in the gain layer becomes positive and continues 
rising because the mass in the gain layer and the energy per 
nucleon grow. The model has passed the critical threshold and 
is on its way to explosion. A closer inspection of the involved 
energies shows that this explosion is powered by neutrino heat¬ 
ing. 

This qualitative difference of the outcome of 2D simu¬ 
lations with 90° and 180° grids is another confirmation of 
the proximity of our 2D simulations, and in particular of 
the 11.2 Mq case, to a success of the convectively supported 
neutrino-driven mechanism. Together with the recent models 
for stars in the 8-10 M© range with 0-Ne-Mg cores, which ex¬ 
plode even in spherical symmetry (Kitaura et al. 2006), our cur¬ 
rent results seem to indicate that the neutrino-heating mecha¬ 
nism is viable at least for stars near the low-mass end of super¬ 
nova progenitors. 

The sensitivity to numerical variations, however, also 
stresses the need to remove some of the shortcomings and limi¬ 
tations of axially symmetric simulations. One must suspect that 
in 3D simulations morphological differences of the structures 
(plumes instead of azimuthal tori), different growth rates of in¬ 
stabilities, or additional degrees of freedom (e.g. triaxial asym¬ 
metries and vortex motion caused by Coriolis forces) might 
lead to sizable quantitative differences which could be crucial 
when collapsing stellar cores are close to the threshold for ex¬ 
plosion. Also the existence of the polar axis of a spherical or 
cylindrical coordinate grid is a potential source of numerical 
uncertainties, because it is a coordinate singularity which is 
impenetrable for approaching fluid flow and thus defines a pre¬ 
ferred grid direction. 

Our results therefore suggest the need to strive for 3D sim¬ 
ulations, preferentially without the disadvantages connected 
with the polar grid axis. The importance of low-mode convec¬ 
tion or low-mode hydrodynamical instabilities as suggested by 
our results implies that such simulations will have to be done 
for the full star and cannot be contrained to a limited wedge. 
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Appendix A: The progenitor modeis 

The properties of the nine progenitor models used in this work 
are summarized in Fig. A. 1 and Table A. 1. In Fig. A. 1 the initial 
pre-collapse data as given by the stellar evolution modelers are 
displayed, in Figs. A.2 and A.3 we compare the cores at a cen¬ 
tral density of pc = 10^^ g cm“^, which was reached by evolv¬ 
ing the models with our ID Boltzmann transport code Vertex. 
Up to this point the infall velocities are still subsonic and neu¬ 
trinos stream off almost unhindered. Compared to the original 
progenitor data in Fig. A.l the electron number has changed 
significantly and strongest for those models which started out 
with the lowest initial densities at the center. The evolution pro¬ 
ceeded nearly adiabatically, the entropy changes are therefore 
small. 

Looking at Fig. A.2, we see that the density structure of 
the inner ~ 1.0 Mq of the iron core is extremely similar in 
all progenitors, correlated with only small differences in the 
electron fraction Fg (Fig- A.3). The central value of Fg varies 
only by 5% between the progenitors. In contrast, the entropy 
per baryon v and infall velocity v exhibit differences of up to 
40%. There is a general trend that |v| and v increase with the 
ZAMS mass, while Fg decreases. 

Outside of the core, larger differences between the progen¬ 
itors exist and are associated with the location of the interfaces 
between layers of different chemical composition and the den¬ 
sity structure in these layers. In particular, the Fe-Si and the 
Si-0 interfaces can have significant influence on the evolution 
of the supernova shock. The enclosed mass at which the inter¬ 
faces are located differs strongly between the progenitors and 
increases non-monotonically with the ZAMS mass, see Table 
A.l. In most cases the composition changes discontinuously 
from the heavy to the lighter nucleus at the interfaces, but also a 
more gradual enrichment of a layer with elements of the neigh¬ 
boring shell is possible. In both cases the progenitor structure 
shows a more or less large entropy jump at the interfaces. For 
larger entropy jumps, which are underlined in Table A.l, also 
large steps occur in the density profile. 

Appendix B: One-dimensional models: more 
details 

With the exception of Model nl3, the collapse times from a 
core central density of pc = 10^^ gcm“^ until bounce dif¬ 
fer only by 10-20% between the models, lying between ~ 20 
and ~ 25 ms (see Fig. B.l and Table B.l). In general, the re¬ 
sults listed in Table B.l and visible in Figs. A.3 and B.1-B.3 
reveal an astonishing degree of convergence of the core evo¬ 
lution for the different progenitors during the phases of col¬ 
lapse, shock formation, prompt shock propagation and break¬ 
out Vg burst (cf. also Liebendorfer et al. 2002). Initially, the 
stellar cores differ in the amount of deleptonization and have 
been evolved differently close to the onset of collapse. Less 
deleptonized cores, however, need a longer time until the col¬ 
lapse becomes dynamical (because their Fg and thus electron 
pressure are higher and photodisintegrations of Fe-group nu¬ 
clei proceed more slowly), thus allowing deleptonization to 
catch up with that of the more evolved progenitors. Despite 


of remaining differences at pg = 10^^ gcm“^ (Figs. A.3 and 
B.2), the central quantities and radial profiles, in particular of 
entropy and Fg, become very similar after neutrino trapping 
sets in at pc ^ 10^^ gcm“^ (Figs. B.2, B.3). This suggests a 
strongly self-regulating character of the hydrodynamics cou¬ 
pled with the neutrino transport, which ensures that the char¬ 
acteristic properties of the homologous core at bounce are al¬ 
most independent of the initial conditions (Liebendorfer et al. 
2002). Only Model s25a28 has such a high value of the entropy 
and such a low value of Fg at the time when pg = 10^^ g cm"^ 
(Figs. A.3 and B.2) that both quantities do not fully converge 
to those of the other models. In this model the higher entropy 
implies that the EoS yields a much larger free proton abun¬ 
dance, which leads to slightly stronger deleptonization and thus 
a somewhat larger entropy (the entropy level is higher by about 
Av ~ 0.07 ke/nucleon) after trapping^. 

As a consequence of the similar collapse evolution, all 
models lose approximately the same amount of energy via neu¬ 
trino emission during collapse, i.e. around 10^^ erg, and be¬ 
cause of the similar core structure and collapse history, the 
shock in all cases is created at an enclosed mass of Mgg ^ 
0.49+0.01 Mq, corresponding to a radius of rgg ^ 10.7+0.1 km 
(Table B.l and Fig. B.3 a). Following Bruenn & Mezzacappa 
(1997) we define the shock creation (sc) time and location by 
the moment and position where the entropy first reaches 3 
per baryon. From now on all times will be normalized to the 
time of core bounce, ^b- 

Also the prompt shock propagation is quite similar. In most 
models, the velocities behind the shock become negative af¬ 
ter 0.9 ms (time ^pse in Table B.l) when the shock encloses a 
mass of Mpsg = 0.78 + 0.01 Mq and has reached a radius of 
Tpsg = 32+1 km (pse stands for “prompt shock ends”). Only in 
case of the progenitors from Limongi et al. (2000) the prompt 
shock pushes out a bit farther. It is interesting to note that in 
all cases the stagnation of the shock happens earlier than the Vg 
burst is released. The corresponding time Le-burst in Table B.l 
is defined as the moment when the Vg luminosity maximum is 
produced at the shock^. 

Although the prompt shock “stalls” in the above defined 
sense, the radial expansion velocity of the shock remains large, 
see Fig. B.4, upper left panel, because of mass of the collapsing 
star being accreted through the shock and accumulating on the 
central, collapsed core. When the shock has passed the neutri- 
nosphere after 4 ms at a radius of 60-70 km (Table B.l), the 
energy and lepton number drain via neutrino emission in the Vg 
burst reduces the thermal and degeneracy pressure of the elec¬ 
trons in the accreted material, so that the shock continues to 
expand more slowly. The strength of this neutrino burst is very 
similar for all models, see Fig. B.4, lower left panel (cf. also 
KachelrieB et al. 2005). 


^ Note that the small differences in Fg and in Fiep after trapping dis¬ 
appear completely, and those in entropy remain slightly larger, when 
the simulations are done with the sophisticated implementation of 
electron captures on nuclei of Langanke et al. (2003). 

^ Technically this moment is determined by taking the time of max¬ 
imum luminosity at 400 km minus the time of flight between the neu- 
trinosphere and this radius. 
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Fig.A.l .The progenitor data for tempera¬ 
ture r, electron fraction Fe, density p, and 
the mass fractions of Si and O, Xsi and Xq, 
respectively, as provided to us by the stellar 
evolution modelers (see Table A.l for ref¬ 
erences). The binding energy and entropy 
were derived with the EoS used in our sim¬ 
ulations. 



r [km] r [km] 


Fig. A.2.Radial profiles of the models dur¬ 
ing core collapse when the central density is 
Pc = 10^^ gcm“^. 


Our models show that the four parameters which character¬ 
ize the subsequent quasi-stationary accretion phase, the mass 
accretion rate through the shock, the mass Mpns and 

radius rpNs of the proto-neutron star (PNS), and the neutrino 
luminosity Ly, are not independent but coupled (Fig. B.4). The 


governing variable is the time-dependent mass accretion rate 
through the shock, which is determined by the progenitor struc¬ 
ture. Since the mass accreted by the shock is further advected 
onto the PNS with a small time delay, the mass of the PNS 
is essentially the time integral of plus an initial value. 
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Fig.A.3 .Profiles as functions of enclosed 
mass of the models during core collapse 
when the central density ispc = 10^^ g cm"^. 
The data are only shown for the stellar mass 
range that was included in the simulations. 
It corresponds to r < 10000 km. 


Starting at ^pb = 4 ms, this initial PNS mass is approximately 
1.0 Mq for all models due to the similar progenitor core struc¬ 
ture. Also the neutrino luminosity depends on The grav¬ 

itational binding energy of the accreted matter must be radiated 
away when accretion proceeds in a stationary way. The total 
neutrino luminosity therefore contains a part from the cool¬ 
ing of the contracting core of the PNS (which loses its bind¬ 
ing energy over a timescale of about 10 s) plus a contribution 
from the settling accretion layer, which is proportional to 5?Msh 
with some lag, because the material has to fall from the shock 
to the cooling layer, where neutrinos are then released over a 
thermal cooling timescale of some 10 ms. Interestingly, also 
the neutrinosphere radii reveal a variability with A large 


(small) accretion rate leads to an approximately stationary so¬ 
lution with more (less) matter piling up on the nascent neutron 
star before it can radiate away its energy in form of neutrinos. 
Thus the PNS obtains a hot, extended (cool, narrow) mantle 
with relatively high (low) densities and optical depths, releas¬ 
ing neutrinos only on a longer (shorter) timescale. The neutri¬ 
nosphere therefore moves to a larger (smaller) radius Vy in case 
of mass accretion proceeding at high (low) rates (Fig. B.4). 

In summary, all the variables which determine the structure 
of the postshock accretion layer depend on the time evolution 
of (9^Msh: A high leads to larger values of Vy and Ly, and 
a faster increase of Mpns- Figure B.4 shows this dependence. 
Looking, e.g., at ^pb = 80 ms we see that dtM^h varies by a 
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Table A.l. List of progenitors used in the simulations. Mpe-si and Msi_o are the enclosed masses at composition interfaces defined by entropy 
jumps. With (e) we denote that a shell interface is connected with a gradual enrichment of the lighter nucleus, i.e. these cases are Fe-FeSi and 
Si-SiO interfaces where the mass fractions of Si or O grows gradually outwards. The underlined numbers indicate an entropy increase of more 
than 1 per nucleon in case of an Fe-Si interface and an entropy increase of more than 2 k^jhy in case of a Si-0 interface. 


Model 

4^zams [^oI 

Reference 

^Fe-Si [^Ol 

Msi-o [Mq] 

Notes 

sib 

>2.26 

Woosley, personal comm. 

1.34 

1.38 

model for Type Ib supernov a 

sll.2 

11.2 

Woosley et al. (2002) 

1.24 

1.30(e) 

from WWW.St ellarevolution.org 

nl3 

13 

Nomoto & Hashimoto (1988) 

1.18 

1.50 


sl5s7b2 

15 

Woosley & Weaver (1995) 

1.28(e) 

1.43(e) 


115 

15 

Limongi et al. (2000) 

1.16(e) 

1.60 


sl5a28 

15 

Heger et al. (2001) 

1.42(e) 

1.81 


s20.0 

20.0 

Woosley et al. (2002) 

1.46 

1.82 

from www.stellarevolution.org 

125 

25 

Limongi et al. (2000) 

1.58 

1.84(e) 


s25a28 

25 

Heger et al. (2001) 

1.62 

1.98(e) 



Table B.l. Characteristic parameters of the ID models for the phases of collapse, prompt shock propagation, and neutrino burst. ?coii is the 
time between the moment when the collapsing core reaches a central density of 10^^ g cm"^ and the moment of shock creation (which is nearly 
identical with the time of core bounce), i.e. the time when the entropy behind the shock first reaches a value of 3 k^fby. The shock creation 
radius and enclosed mass Msc are defined by the radial position where this happens. The energy loss via neutrinos during the collapse 
phase is evaluated by integrating the total neutrino luminosity (for an observer at rest at r = 400 km) over time from the moment when the core 
reaches pc = 10^^ gcm“^ until the moment when the dip in the Ve luminosity is produced around 2.5 ms after shock formation. We call the 
time when the velocities behind the shock drop below 10^ cms"^ the end of the prompt shock propagation phase. This time, tpse, is measured 
relative to the moment of shock creation. At the end of the prompt shock propagation phase, the shock is at radius rpse and its enclosed mass 
is Mpse- The time of the Ve burst, Le-bursts is defined as the post-bounce time when the Ve luminosity maximum is produced at the shock, which 
then is located at the radius -burst and mass -burst- Finally, the energy emitted during the prompt Ve burst is defined as the time-integrated 
Ve luminosity for the FWHM of the burst, L'^ursf evaluated at 400 km for an observer at rest. 


Model 

Loll 

[ms] 

''sc 

[km] 

[Mq] 

77V,loss 

^coll 

[10^^ erg] 

^pse 

[ms] 

^pse 

[km] 

^pse 

[Mq] 

^Ve-burst 

[ms] 

^ve -burst 

[km] 

^Ve-burst 

[Mq] 

^burst 

[lO^i erg] 

sib 

23.9 

10.7 

0.49 

1.00 

0.87 

32 

0.78 

3.8 

64 

1.00 

1.48 

sll.2 

25.2 

10.7 

0.50 

1.00 

0.87 

32 

0.78 

3.7 

63 

1.00 

1.38 

nl3 

28.9 

10.7 

0.50 

0.96 

0.88 

32 

0.78 

3.7 

61 

0.98 

1.29 

sl5s7b2 ^ 

23.9 

10.7 

0.49 

1.01 

0.91 

32 

0.78 

3.8 

64 

1.00 

1.47 

115 

20.9 

10.6 

0.50 

1.04 

1.03 

35 

0.82 

4.1 

68 

1.04 

1.74 

sl5a28 

21.6 

10.7 

0.49 

1.03 

0.95 

33 

0.79 

4.0 

66 

1.02 

1.60 

s20.0 

22.7 

10.6 

0.49 

1.03 

0.91 

32 

0.79 

3.8 

64 

1.01 

1.50 

125 

20.0 

10.7 

0.49 

1.06 

1.15 

37 

0.83 

4.1 

68 

1.03 

1.82 

s25a28 

19.8 

10.6 

0.48 

1.06 

1.00 

34 

0.79 

4.1 

68 

1.03 

1.93 


^ This model is identical with Model sl5Gio_ld.b in Buras et al. (2006) except for a slightly different, improved implementation of the 
relativistic corrections to the gravitational potential. This difference causes only small changes in the neutrino luminosities of a few percent, 

but otherwise has no visible consequences. 



0.0 -0.5 -1.0 -1.5 -2.0 -2.5 -3.0 -3.5 

log(teb-t) [s] 


Fig. B.l. Central density versus time remaining until core bounce for 
all ID models. 


factor of five between the models, being 0.8 and 1.3 Mq 
in case of Models nl3 and si 1.2, respectively, at the low mass 
end of considered progenitors, and 3-4 Mq/s for Models 125 
and s25a28 at the high end of the progenitor mass spectrum. 
Consistently, ry at ^pb = 80 ms is only 60 km for the low-mass 
progenitors and 25% higher, i.e. about 80 km, for the high-mass 
models. Also the luminosities differ by a factor of two, and the 
PNS masses show differences of order 25%. 

The shock radius follows the time evolution of the neu¬ 
tron star radius qualitatively, but with some time delay. 
Interestingly, during phases of slowly changing the influ¬ 
ence of different values of dtM^h on the shock position is rather 
modest. The shock radii are nearly identical for all progenitors 
until ~ 25 ms after bounce and differ at most by about 10% until 
about 80 ms post-bounce. During this phase, the significantly 
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Fig.B.2. Central entropy (panel a) and central electron and lepton 
fraction (panel b) versus central density during core collapse for all 
ID models. The vertical dotted line marks the time of comparison at a 
central density 10^^ gcm"^ in Figs. A.2 and A.3. 

larger differences in Mpns, and Ly mentioned above 

seem to partly compensate each other in their influence on the 
shock radius. For higher values of the shock radius is 

smaller because of the higher ram pressure. A larger Mpns in¬ 
creases the gravitational pull and thus also lowers But on 
the other hand a higher Ly increases the heating and thus leads 
to an expansion of the heating region, which is supported from 
below by a larger PNS (because Vy is larger). 

Of course, the counteracting effects do not compensate each 
other perfectly. As can be seen at later times, the shock radii 
differ more strongly when the PNSs begin to show a wider 
spread in mass, although with nearly the same radius, and in 
particular when the mass accretion rates show differences of a 
factor of about 10 instead of the factor of 2 at early times. The 
analytic study of Arcones Segovia (2003) and Arcones & Janka 
(in preparation) reveals that the steady-state accretion shock ra¬ 
dius is a sensitive function of the PNS radius (r^h ^ ^y^ '^), 
and decreases less strongly with higher mass accretion rate and 
proto-neutron star mass (see also Fryer et al. 1996). 

Note that a sudden decrease of as it happens when 
an infalling composition shell interfaces reaches the shock (see 
Sect. 2, also nicely demonstrates the dependence of ry, Ly, and 
MpNs on this variable: Immediately after such a drop, Ly starts 
decreasing, the decrease of ry becomes steeper, and the increase 
of MpNs slows down (Fig. B.4). Also, the total lepton num¬ 
ber and energy loss via neutrino emission, see Fig. B.5, de¬ 
pends strongly on dtM^h and thus on the progenitor structure 



b ^ [^^®] 


Fig. B.3. Entropy (panel a) and electron and lepton fraction (panel b) 
versus enclosed mass at the time of shock formation. 

during the post-bounce accretion phase. These facts could be 
used to infer the progenitor structure from a supernova neu¬ 
trino signal which is measured with good time resolution (cf. 
also Liebendorfer et al. 2003). 

Appendix C: Linear anaiysis of PNS structurai 
changes 

A linear analysis is performed to determine the sensitivity of 
the PNS structure, i.e. contraction or expansion, to changes or 
redistribution of lepton number and energy in some layer. This 
analysis is then applied to the main effects of convective ac¬ 
tivity in the nascent neutron star, (1) the transport of lepton 
number and energy within the convective layer, and (2) the dif¬ 
ferences in the loss of lepton number and energy by neutrino 
emission in 2D simulations compared to ID models. 

We demonstrate the analysis in case of lepton number vari¬ 
ations, holding the specific internal energy fixed, but variations 
of the internal energy can be investigated in the same way. 
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Fig.B.4. Post-bounce evolution of the models in spherical symmetry. The left panels show phases around the core bounce with the last stage 
of infall it < 0), prompt shock propagation, neutrino burst, and early accretion phase of the shock. The right panels depict the post-bounce 
accretion phase of the stalled shock. From top to bottom, the panels display the shock radius, the mass accretion rate through the shock, the mass 
enclosed by the shock (left) or PNS baryonic mass (right), the PNS radius, and the electron neutrino (left) or total neutrino (right) luminosity at 
r - 400 km for an observer at rest. The PNS mass and radius are defined by the corresponding values at the Ve sphere. 


Assuming Newtonian gravity, which is sufficient here because with G being the Newtonian gravitational constant. Using the 
we are interested in qualitative, not quantitative results, the hy- relation dM = 47ipr^dr , we get 
drostatic structure of the PNS is determined by 


dp(r) GM(r)p 

pi 


dp{M) GM 

(c-i) -^ = - 


(C.2) 


dr 


Anr^iM) ' 
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Fig.B.5. Integrated lepton number (top) and energy loss of the ID 
models versus time. 


An infinitesimal change SY(M) of lepton number, e.g. by re¬ 
moving one electron from a layer with infinitesimal mass AM 
and baryon number AAby at mass coordinate Mp (we here and 
in the following drop the subscript of Fiep), i.e. 

6Y(M) = -1 • A;V-' = - (AM/mby)”' (C.3) 

if Mp < M < Mp -I- AM, and SY(M) = 0 otherwise, will lead to 
a variation of the pressure according to 


«M)=| 


mmM)+ & 

I uf) 


(M) Sp(M ), 


(C.4) 


where Sp(M) is the structural change induced by the change of 
lepton number. The new hydrostatic solution is 

djp^Sp) ^ GM 
dM 47r(r -l- SrY ’ 

where Sr(M) is the movement of mass shell M. 

Expanding the RHS to first order, inserting Eq. (C.4), and 
subtracting Eq. (C.2) leads to 

GMSr 


d 

dp 

dp 

•SY+ — 

‘Sp 

dM 


P dp 

Y 


r^5 


(C.6) 


Now setting the RHS of Eq. (C.6) to zero for the moment, the 
solution reads 


where we have expanded to first order and made use of 
Eqs. (C.7) and (C.3) and of the relation AM = 47Tr^p(M^)Sr. 
The change in mass SM is independent of R and means expan¬ 
sion, if negative (i.e. less mass within a given radius R). 

A relation between changes of mass and radius can be ob¬ 
tained by using M = 47rpr'^dr' = 47r(p + Sp)/^d/ . 

Expanding to first order, we thus get a relation between Sr and 
Sp: 

1 

Sr = - - I 47TSpr'^dr' . (C.9) 

47rpr2 Jo 

Note that the radius displacement SR{M) depends on R, see 
Eq. (C.9), and is therefore less suitable for our analysis. 

In order to take into account the RHS of Eq. (C.6) we apply 
an iterative procedure: The solution Sri(M) of iteration step i 
is inserted on the RHS of Eq. (C.6), which is then numerically 
integrated over M to find the solution Spi+i{M). In the first it¬ 
eration step the solution for Sp{M) from Eq. (C.7) is used to 
calculate Sri(M) from Eq. (C.9). We find that for typical PNS 
profiles Eq. (C.7) provides a very good approximation already; 
the next iteration step changes SM by 10“^ at most. 

In Eig. C.I we present the change of SM when one elec¬ 
tron or an energy of 10 MeV is removed from the PNS at dif¬ 
ferent radii (or densities, corresponding to different enclosed 
masses Mp) for two representative times in the evolution of 
Model sI5_32. Note that adding one electron or 10 MeV would 
have the opposite effect. We see that removing this amount of 
energy always leads to contraction, while the extraction of one 
electron always causes expansion of the PNS. The latter finding 
can be understood from the fact that the electron is taken away 
but its degeneracy energy is assumed to remain in the PNS as 
thermal energy. 

When discussing the effects of the convective transport of 
energy and lepton number on the PNS structure, it is interesting 
to compare the values of SM in the drain region, where ^int and 
Fiep drop, and in the dump region, where Cint and Fiep increase 
due to convection. For example, at r = 73 ms, the removal of 
10 MeV from the drain region will lead to an increase of the 
mass enclosed by some large radius R by (0.7-1.3) x lO'^'^g 
(Fig. C.I), while adding these 10 MeV in the dump region will 
reduce SM(R) by (0.9-2.1) x lO'^'^g (Fig. C.I). This means a 
net decrease of SM(R), and thus an expansion of the PNS. This 
behaviour is also seen at later times and in a similar way for the 
transport of Fiep. 

Appendix D: Mixing scheme for reproducing PNS 
convection 


dp 


(M). 


dp 

dp 


(M) 


6Y{M). 


(C.7) 


In terms of the “mass shift” 6M the resulting expan¬ 
sion/contraction of the PNS is now evaluated at a radius R > 
r(Mp): 

r-R 


6M{R) 


T 


47r^p(M(r'))r'Mr' 




dp 

dp 


(Mp) 


^by 

' p(Mp) ’ 


(C.8) 


In order to test the consequences of PNS convection on the PNS 
and supernova evolution in ID simulations, we have developed 
a simple numerical “mixing algorithm”, which reproduces the 
energy and lepton number transport found in 2D models to a 
high degree. 

In this algorithm we assume that any convective activity 
in the PNS leads instantaneously to a redistribution of energy 
and lepton number in the unstable region so that the convec- 
tively unstable gradients disappear, i.e. Cql = 0 is estab¬ 
lished. Setting y^diff = 1 and assuming d^Fiep^ /dr = dFiep /dr 
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Fig.C.1 .Results of a linear analysis of the 
structural changes of the PNS in response 
to the removal of lepton number or energy, 
respectively, at different densities and thus 
enclosed masses Mp. The lower lines cor¬ 
respond to the expansion when one elec¬ 
tron is removed, while the upper lines depict 
the contraction when an energy of 10 MeV 
is extracted. The lines correspond to post¬ 
bounce times of 73 ms (solid) and 240 ms 
(dashed). The labels “dump” and “drain” 
denote the density windows between the 
vertical, dotted lines, which roughly indi¬ 
cate where PNS convection deposits or ex¬ 
tracts, respectively, lepton number and en¬ 
ergy. 



a [Mo] 

tpb= 46.5ms 



r [km] r [km] 

Fig.D.l. Panel a depicts the Fiep profiles for three times around the 
onset of PNS convection, comparing a ID model (sl5s7b2), a 2D 
model (sl5_32, for which angular-averaged values are plotted), and 
a ID model with PNS mixing scheme (sl5_mix). Panel b shows Plep, 
the total entropy the radial velocity v^, and the Brunt-Vaisala- 

frequency 6 c>>bv (clockwise from top left to bottom left) in the PNS 
convective region for Model sl5_32 shortly after the onset of PNS 
convection. 


in Eq. (7), we find that the instability is mainly driven by a 
negative gradient of the entropy (including neutrino entropy). 
Cql = 0 then corresponds to a flat entropy profile, consistent 
with what is observed in the PNS convective region in 2D sim¬ 
ulations, see Figs. 9 and 10. The numerical mixing scheme re¬ 
distributes the energy in an unstable layer in a conservative way 
such that the entropy profile develops a plateau. 

As a second constraint we introduce the empirically derived 
relation 


T/ new 
Tlep 




1 a 


^new _ ^old 
^int ^int 


old 

int 


(D.l) 


where “old” denotes the lepton number Fiep and internal energy 
Cint before the application of the mixing scheme in a time step 
of the simulation, “new” denotes the variables after such an ap¬ 
plication, and A is a global normalization factor which ensures 
that lepton number is globally conserved. 

The mixing algorithm is executed in an operator splitting 
step after each hydrodynamic step. The algorithm detects re¬ 
gions with a negative entropy (s-\-Sy) gradient inside the PNS. In 
these regions energy and lepton number are then redistributed 
in such a way that Eq. (D.l), d(s -f Sy)/dr = 0, and global 
energy and lepton number conservation are fulfilled (GR grav¬ 
itational effects are ignored). 

We find empirically that for a = 2 the evolution of the pro¬ 
files of Fiep, ‘S' -f and ^int inside the PNS for simulations with 
the mixing scheme reproduce very well the results of the 2D 
model, see Fig. 10. Near the boundaries of the convectively 
unstable layer small differences can appear. The over- and un¬ 
dershooting of dynamically moving matter into convectively 
stable layers can, of course, not be properly accounted for by 
this simple scheme. Another deficiency of the scheme is best 
visible shortly after the onset of convection around 50 ms af¬ 
ter bounce, when the Brunt-Vaisala frequency can reach val¬ 
ues of up to 2 ms“^ (Fig. 9 and Fig. D.l, panel b), and the 
convective layer has not yet fully developed and the transient 
starting phase with growing perturbations is not yet over. In 
this nonstationary, early phase of PNS convection the results 
with the mixing scheme cannot reproduce the Fiep-profiles of 
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M [Me] 

Fig.E.l. Maximum fluctuations of the density within the core of 
1.5 Mq at four times during the collapse of Model sl5_64_p (top). The 
initially imposed random perturbations are damped strongly during 
the subsonic collapse in the inner part of the Fe core. The correspond¬ 
ing density and velocity profiles (averaged in latitude) are shown in 
the middle and bottom panels, respectively. The short solid lines cut¬ 
ting the profiles in the bottom plot indicate the position of the sonic 
point. 

the 2D simulations very well (Fig. D.l, panel a); lepton num¬ 
ber is transported more efficiently in the 2D models. 

Appendix E: Perturbation growth during coiiapse 

The majority of our models were computed through the phases 
of core collapse, bounce, shock formation, and shock stagna¬ 
tion in spherical symmetry, and then continued in 2D with ran¬ 
dom perturbations added during the mapping from the ID to 
the 2D grid. In contrast, we followed the evolution of Models 
sl5_64_p and sl5_64_r in two dimensions also during core col¬ 
lapse, shock formation, and early shock propagation. In this 
case the spherically symmetric density structure of the progen¬ 
itor core was perturbed in each cell of the computational grid 
in a random way, in case of Model sl5_64_p with an amplitude 
of 2%, in Model sl5_64_r with 1% amplitude (Table 1). We 
will discuss here briefly the evolution of the perturbations dur¬ 
ing the core collapse phase, concentrating on Model sl5_64_p. 



Fig. E.2. Standard deviations of the density fluctuations in lateral di¬ 
rection, cTp, as defined in Buras et al. (2006, Eq. 24) for different mass 
shells versus the (time-dependent) radius km of these mass shells. The 
dashed lines represent fits according to cTp oc the thick black bars 
mark the phases of the collapse when the infall velocity of a mass shell 
is more than 60% of the local sound speed, and the white bars mark 
the phases when the local infall velocity is supersonic. 


Since the initial spin is too low, rotation hardly affects the pre¬ 
bounce evolution of Model sl5_64_r, and the latter model does 
not show any important rotation-induced differences with re¬ 
spect to the perturbation growth. 

Because of the lack of a force continuously driving their 
creation, the initial perturbations are rapidly damped during 
the still subsonic early phase of the collapse. This damping 
is stronger in the central regions of the core where the sound 
speed has the highest values and the timescales of sonic cou¬ 
pling are shortest (Fig. E.l, top panel). When supersonic col¬ 
lapse velocities are reached, which is the case only about 10 ms 
prior to bounce (Fig. E.l, lower panel), even in the outer 
regions of the core the maximum density fluctuations have 
dropped from initially 2% to less than 1%. Only at the outer 
edge and outside of the iron core (enclosed masses of more 
than 1.2 Mq) the densities and sound speed are so low that the 
maxima of the density fluctuations are still near their initial 
values. Once the collapse velocities reach the sound speed and 
become supersonic, the standard deviations of the density fluc¬ 
tuations begin to grow again slowly like the inverse square root 
of the decreasing radius (cTp oc Fig. E.2), following the 

prediction by Lai & Goldreich (2000). The peak fluctuations, 
however, have no time to rise again significantly until bounce 
(Fig.E.l). 

When running the simulations continuously in 2D instead 
of mapping ID post-bounce models on a 2D grid some mil¬ 
liseconds after bounce, we detect a slightly earlier onset of 
convection, which is seeded by the density and velocity fluc¬ 
tuations in the convectively unstable layers that develop below 
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the neutrino sphere and in the neutrino-heating layer behind the 
stalled shock. Despite the earlier start of convection behind the 
shock, no significant differences of the evolution could be dis¬ 
covered between Models sl5_32 and sl5_64_p (Figs. 11-17). 
We therefore have included the results of the latter model in 
the relevant plots in Sect. 3, but have not specifically discussed 
them in the text. 



